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ABSTRACT 


Watershed  modelling  of  gaged  streams  in  the  Chemung 

and  Upper  Susquehanna  River  Basins  was  attempted  using  the 
U.S. Array  Corps  of  Engineers'  HEC-1,  Flood  Hydrograph  Pack¬ 
age.  The  Clark  method  of  unit  hydrograph  theory  and  the 
initial  and  constant  loss  rate  approach  were  used  in  this 
analysis.  Parameters  were  derived  from  a  total  of  52 
storm  events  for  nine  watersheds.  Regression  equations 
were  derived  relating  the  Clark  unit  hydrograph  parameter 
to  physical  basin  characteristics.  Regional  constant  loss 
rates  were  established. 
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CHAPTER  ONE 
INTRODUCTION 


Floods  are  headline  grabbers.  Every  year  floods  show 
us  how  destructively  nature  can  unleash  its  power.  They 

receive  extensive  news  coverage  which  describes  the 
resulting  damage  and  loss  of  life.  Floods  disrupt  lives 

and  productive  endeavors.  A  great  deal  of  effort  has 
gone  into  studies  of  past  floods  and  estimates  of  future 
floods  in  order  to  develop  plans  for  flood  damage  re¬ 
duction  and  prevention.  Richards  (1955)  wrote,  "The  flood 
problem  is  one  which  affects  many  branches  of  engineering 
other  than  those  dealing  primary  with  the  control  of 
water."  He  also  wrote  that  the  flood  estimation  problem 
is  a  peculiarly  difficult  and  complicated  one.  Although 
the  natural  laws  governing  floods  are  both  recognizable 
and  generally  appreciated,  the  difficulty  lies  in  the 
application. 

The  flood  factors  which  are  normally  considered  the 
most  important  are  as  follows i  rain  intensity,  rainfall 
distribution,  rainfall  duration,  drainage  area,  basin 
shape,  basin  slope,  land  cover,  and,  finally,  the  initial 
state  of  wetness  in  the  soil  (Richards,  1955,  P*  8)* 

Hoyt  and  Langbein  (1955)  recognized  the  need  for  a 
national  policy  addressing  the  flood  problem.  They  hypothe¬ 
sized  that  some  future  Congress  or  Congresses  would  recog- 


1 


2 


nize  that  "a  sound  approach  to  floods  must  be  broader 
than  flood  protection  alone,...".  In  their  ten  point 
proposal  number  three  is  quoted,  "Flood-forecasting  ser¬ 
vices  should  be  recognized  as  an  essential  part  of  flood 
management  and  be  appropriately  supported."  Indeed,  around 

1955*  the  U.S.  Weather  Bureau  (now  the  National  Weather  Ser¬ 
vice,  NOAA)  received  appropriations  for  its  flood- forecast¬ 
ing  service  averaging  about  $900, 000/year .  The  savings 
accredited  to  this  service  were  roughly  estimated  at  about 
$27,000,000,  implying  a  cost-benefit  ratio  of  1  to  30. 

Flood- forecasting  is  certainly  not  a  modem  notion. 
Diodorus  Siculus  visited  Egypt  during  the  1st  Century  B.C. 
(Biswas,  1967,  p.125).  On  the  Nile  River,  the  stage  of  the 
Nile  were  maintained  by  using  nilometers  (gaging  stations). 
These  structures  were  used  to  observe  the  river  stage  and 
compare  the  current  height  with  that  of  previous  high 
water  marks.  Little  evidence  of  actual  technical  co-opera¬ 
tion  between  the  temples  charged  with  running  the  nilometers 
exists  but  according  to  Diodurus,  flood  warnings  were 
sounded  to  the  population  from  the  nilometer  at  Memphis  (the 
seat  of  the  government)  in  case  of  emergency.  The  ancients' 
warning  system  worked  as  follows;  with  the  approach  of 
the  flood  season,  the  stage  of  the  Nile  was  carefully 
watched  and  compared  with  the  markings  of  the  previous 
year.  Swift  rowers  were  sent  from  the  furtherest  upstream 
gaging  stations,  one  after  another,  to  report  the  latest 
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level  at  the  capital.  These  extremely  good  rowers,  rowing 
with  the  current,  were  able  to  outpace  the  approaching 
peak  flood  and  give  sufficient  advance  warning  to  the 
people  along  the  river  of  the  forthcoming  catastrophe. 

No  one  would  seriously  consider  instituting  a  flood 
prediction  system  today  based  on  a  fleet  of  rowboats.  How¬ 
ever,  there  is  a  need  to  forecast  floods  in  a  "real¬ 
time"  setting  in  order  for  warnings  to  be  broadcasted, 
people  and  property  to  be  evacuated,  and  emergency  re¬ 
sources  to  be  mobilized  in  the  area  of  the  expected  threat. 

But  how  does  one  go  about  forecasting  a  flood?  This 
thesis  will  present  only  one  method.  That  method  is  the 
unit  hydrograph  approach  for  ungaged  watersheds  which 
was  suggested  by  the  U .S .  Army  Corps  of  Engineers  in  its 
manuals  and  which  will  be  presented  later  in  Chapter  Three . 
The  selection  of  this  method  was  based  in  part  on  the 
observation  of  Dooge  (1959)  who  wrote,  "The  unit  hydrograph 
approach  to  stream  flow  has  developed  into  one  of  the 
most  powerful  tools  of  applied  hydrology." 

GOAL  AND  SCOPE  OF  THESIS 

The  goal  of  this  thesis  is  to  use  the  U.S.  Army 
Corps  of  Engineer's  HEC-1  computer  program  (discussed  in 
Chapter  Three)  in  developing,  (a)  Regression  relationships 
between  the  unit  hydrograph  parameters  (based  upon  the 
Clark  Method)  and  some  of  the  more  important  physical 


basin  characteristics;  and,  (b)  Regional  relationships  to 
determine  the  initial  loss  and  constant  loss  rate  parameters 
in  .the  area  of  interest. 

The  scope  of  this  thesis  is  to  accomplish  the  goal 
by  using  historical  storm  and  rainfall  data  from  several 
gaged  streams  in  the  area  of  interest  in  solving  the  inverse 
problem.  The  inverse  problem  is  defined  as  a  problem  in 
which  the  input  (rainfall)  and  the  output  (stream  runoff) 
are  known  and  are  used  in  determining  the  transfer  functions, 
i.e.  the  Clark  unit  hydrograph  and  loss  rate  parameters. 

The  area  of  interest  is  that  of  the  Chemung  and  Upper  Sus¬ 
quehanna  River  basins  in  the  state  of  New  York  and  a  small 
portion  of  Pennsylvania.  Nine  gaged  stream  basins  in  this 
area  were  selected  for  use  in  calibrating  the  computer 
model  HEC-1 .  These  basins  which  are  unregulated  were  chosen 
on  the  basis  of  their  hydrological  similarity.  The  ungaged 
watersheds  in  the  area  around  Coming  and  Elmira,  New  York, 
are  the  main  targets  of  this  modelling  effort. 

OVERVIEW 

This  work  is  organized  into  six  chapters.  Chapter 
One  is  the  introduction  in  which  the  consequences  of  floods, 
the  need  for  flood  prediction,  and  the  goal  and  scope  of 
this  thesis  were  presented.  Chapter  Two  is  a  literature 
review  primarily  concerned  with  the  origins  of  the  HEC-1 


approach  selected  for  this  study.  Chapter  Three  describes 
the  specific  method  employed  to  obtain  the  results. 

Chapter  Pour  presents  the  data  and  calculations  performed 
in  the  application  of  the  method  in  the  area  of  interest. 
Chapter  Five  is  a  discussion  of  the  results.  Chapter  Six 
presents  conclusions  reached  by  this  author  and  other 
authors.  Several  appendices  are  provided  to  keep  the  basic 

chapters  as  free  as  possible  from  "data  overload”  and 
are  identified  in  the  basic  work. 


CHAPTER  TWO 


LITERATURE  REVIEW 

PURPOSE  AND  SCOPE 

The  purpose  of  this  chapter  is  to  trace  the  development 
of  the  unit  hydrograph  theory  and  its  application  in  flood 
prediction  on  ungaged  basins.  The  scope  of  this  chapter  is 
confined  to  a  brief  discussion  of  the  origins  of  the  three 
components  of  the  HEC-1  model  with  a  limited  presentation 
of  other  methods . 


LOSSES 

The  term  precipitation  loss  is  defined  in  this  thesis 
as  the  amount  of  precipitation  not  available  during  a  storm 
for  overland  flow.  There  are  basically  two  types  of  precip¬ 
itation  losses.  One  is  called  interception.  It  represents 
the  amount  of  rainfall  which  is  held  by  the  trees,  grass 
and  local  depressions  such  as  potholes  and  cracks.  Any¬ 
thing  that  traps  and  does  not  allow  the  rainfall  to  move 
about  freely  is  classified  under  this  type  of  loss.  The 
other  type  of  precipitation  loss  is  infiltration.  It  rep¬ 
resents  the  amount  of  rainfall  absorbed  into  the  earth. 

In  computing  rainfall  excess,  the  precipitation  losses 
are  subtracted  from  the  total  storm  rainfall.  Sokolov  et 
al.  (1976,  p.175)  listed  five  methods  to  determine  rain¬ 
fall  losses.  They  are  mentioned  below  1 
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1.  Constant  runoff  coefficient,  which  is  the  ratio  of 
surface-runoff  depth  to  total  rainfall  depth. 

2.  Constant  infiltration  rate  during  a  storm  known  as  the 
-index  in  the  U.S. 

3.  Initial  rainfall  loss  followed  by  a  constant  infiltra¬ 
tion  rate.  In  this  method  all  rainfall  is  considered  to  be 
lost  until  an  initial  rainfall- loss  demand  has  been  satisfied. 
Thereafter,  rainfall  is  lost  at  a  constant  rate. 

4.  Regional  curves  of  infiltration  capacity,  that  are 
variable  with  time,  and  are  based  on  soil  type,  land  use, 
and  an  antecedent-rainfall  index. 

5*  A  variable  runoff  coefficient  CTt  method,  which  varies 
with  time,  even  during  storm  periods,  in  response  to  the 
variation  of  a  basin  moisture  index. 

Any  one  of  these  methods  may  be  included  in  a  basin 
model.  HEC-1  uses  the  methods  described  in  3.  and  4. 
above  and  will  be  discussed  further  in  Chapter  Three. 

UNIT  HYDROGRAPH 

ORIGINAL  UNIT  HYDROGRAPH  THEORY 

Sherman  (1932,  p.  501)  first  proposed  the  unit-graph 
(as  it  was  originally  called)  method  as  a  way  to  predict 
flood-peak  discharges  and  discharge  hydrographs  from  rain¬ 
fall  events.  Sherman  defined  the  unit-graph  as  follows: 

If  a  given  one-day  rainfall  produces  a  1-inch 
depth  of  rainfall  over  the  given  drainage  area,  the 
toyarograph  showing  the  rates  at  which  the  runoff  occur¬ 
red  can  be  considered  a  unit-graph  for  the  watershed. 
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The  unit  hydrograph  (unit-graph)  theory  makes  two 
assumptions!  1)  Runoff  response  to  rainfall  is  time  invar¬ 
iant.  2)  Runoff  response  to  rainfall  is  linear;  this 
implies  the  principle  of  superposition.  These  assumptions 
are  not  rigorous  and  there  is  some  flexibility. 

By  means  of  the  second  assumption,  Sherman  observed 
the  possibility  of  computing  the  runoff  history  correspond¬ 
ing  to  a  rainfall  of  any  duration  or  degree  of  intensity  for 
the  same  watershed  from  a  single  observed  hydrograph.  In 
applying  his  method,  he  noted,  that  the  ordinate  and  time 
intervals  of  unit  hydrographs  for  two  similar  watersheds  of 
different  sizes  are  proportional  to  the  square  root  of  the 
watershed  areas,  provided  the  difference  in  size  is  not 
great.  He  also  observed  that  consistent  results  in  the  data 
were  based  upon  the  observations  being  confined  to  a  single 
area  or  to  closely  similar  areas  and  being  segregated  accod- 
ing  to  the  seasons. 

Unit  hydrograph  methods  are  usually  based  on  the  ability 
to  Separate  the  surface-runoff  hydrograph  from  the  hydrograph 
of  the  total  rainfall.  The  total  hydrograph  is  often  assumed 
to  consist  of  three  types  of  flow  defined  below  by  Dooge 

(1959*  P*  24l)i  1)  Surface  runoff,  or  the  water  reaching 

surface  channels  by  the  overland  route;  2)  Interflow,  or 
the  portion  of  infiltrated  water  that  passes  through  the 
shallower  horizons  of  the  soil  to  reach  defined  stream 
channels  within  a  relatively  short  time,  without  first 
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reaching  the  main  ground-water  table;  and,  3)  Base  flow, 
or  the  water  contributed  as  ground  water  outflow  from  an 
aquifer.  It  should  be  noted  that  interflow  is  usually  com¬ 
bined  with  surface  runoff  and  collectively  called  surface 
runoff. 

LATER  UNIT  HYDROGRAPH  DEVELOPMENTS 

In  193^.  1936  and  1937  Zoch  (Dooge,  1959.  p.24l)  pub¬ 
lished  several  papers  on  the  unit  hydrograph  method  based 
on  the  assumption  that  at  any  time  the  rate  of  discharge  is 
proportional  to  the  amount  of  rainfall  remaining  with  the 
soil  at  that  time.  The  S-hydrograph  method  developed  by 
Morgan  and  Hullinhors  (1939)  employed  a  unit  hydrograph  of 
a  certain  duration  to  form  an  S-hydrograph  resulting  from 
a  continuous  applied  rainfall  to  construct  unit  hydrographs 
for  other  than  the  original  duration. 

Snyder  (1938.  p.^47)  and  Taylor  and  Schwartz  (1952,  p. 
235)  applied  unit  hydrograph  theory  to  basins  by  relating 
features  of  the  unit  hydrograph  to  watershed  characteristics 
by  strictly  statistical  considerations.  Their  attempts  led 
to  the  "synthetic  unit  hydrograph  method." 

In  contrast  to  this  statistical  approach,  several 
workers  attempted  to  give  a  more  conceptual  basis  to  the  unit 
hydrograph.  Clark  (19^5.  P*  1^19)  developed  the  instantaneous 
unit  hydrograph  (IUH)  method  by  suggesting  the  unit  hydro- 
graph  for  instantaneous  rainfall  could  be  derived  by  routing 
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the  time-area-concentration  curve  through  a  single  element 
of  linear  reservoir  storage.  O’Kelly  (1955»  p«  3^5)  of  the 
Irish  Office  of  Public  Works  suggested  replacement  of  the 
time-area-concentration  curve  by  an  isosceles  triangle. 

Nash  (1957t  p.  114)  noted  that  the  instantaneous  unit  hydro¬ 
graph  could  be  derived  by  routing  the  instantaneous  rain¬ 
fall  through  a  series  of  successive  linear  reserviors  of 
equal  delay  time.  Dooge  (1959t  p.  24l)  attempted  to  remove 
many  of  the  subjective  elements  from  unit  hydrograph  analysis. 
The  number  of  IUH  investigations  goes  on  and  on.  Chow  (1964, 
p.  14-25)  compiled  an  extensive  list  of  past  accomplish¬ 
ments  in  his  handbook.  He  concluded  that  the  use  of  the 
IUH  rather  than  other  unit  hydrograph  methods  is  better 
suited  for  investigations  on  the  rainfall  and  runoff  rela¬ 
tionship  in  drainage  basins. 


CLARK’S  METHOD  OP  UNIT  HYDROGRAPH  ANALYSIS 

Clark  (1945*  p.  1422)  was  the  first  to  adapt  the  idea 
of  the  IUH  in  hydrograph  analysis  to  the  unit  hydrograph 
theory.  In  his  original  work  he  wrote  1 

"The  range  of  possible  unit-graph  determinations 
can  be  reduced  by  correlating  recognized  discharge 
concepts  with  the  physical  limitations  of  valley  stor¬ 
age  and  with  the  time  elements  which  necessarily  result 
from  storage  and  discharge  relations .  The  correlation 
can  be  utilized  to  develop,  from  time-area-concentra¬ 
tion  curves  of  specific  drainage  areas,  unit  hydrographs 
with  a  small  range  of  determination  variability,  inde¬ 
pendent  of  assumptions  regarding  runoff  distribution 
in  the  flood-producing  storm  and  reflecting  influences 
of  drainage  area  shape  and  stream  pattern." 
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He  showed  the  relationship  between  the  unit  hydro - 
graph  and  the  methods  of  flood  routing  (means  of  modifying 
a  hydrograph  by  the  effects  of  valley  storage).  Clark 
stated  that  one  advantage  of  using  an  "instantaneous  hydro¬ 
graph  is  that  it  can  be  derived  from  the  fundamental  char¬ 
acteristics  of  the  basin  and  then  used  with  any  length  of 
unit  period  to  determine  the  unit  hydrograph,  Other 

advantages  stated  by  Clark  which  make  his  method  appealing 
to  "real-time"  flood  forecasting  arei  1)  The  procedure  is 
definable.  Identical  hydrographs  will  be  obtained  by  dif¬ 
ferent  people.  2)  Determination  is  independent  of  any 
knowledge  of  runoff  distribution,  except  the  time  of  ending. 
3)  The  unit  hydrograph  quantities  are  instantaneous  rates 
of  discharge  at  the  time  specified. 

FLOOD  ESTIMATION i  UNIT  HYDROGRAPH 

Richards  (19 55»  p.10)  presented  the  principal  factors 
affecting  floods.  A  number  of  flood  formulas  is  discussed 
and  a  brief  outline  of  the  flood  frequency,  probability, 
and  unit  hydrograph  methods  used  in  the  U.S.  is  included. 
Sokolov  (1976,  p.182)  states  that  flood  flow  computations 
for  unit  hydrograph  derivation  are  accomplished  using 
three  general  methods  1  1)  Analysis  of  rainfall — runoff 
records  for  isolated  unit  storms,  2)  Analysis  of  rainfall — 
runoff  records  for  major  storms,  and  3)  Computation  of 
synthetic  unit  hydrographs  from  a)  Direct  analogy  with 
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basins  of  similar  characteristics,  or  b)  Indirect  analogy 
with  a  large  number  of  other  basins  through  the  appli¬ 
cation  of  empirical  relations. 

BASIN  CHARACTERISTICS  USED  IN  HYDROGRAPH  ANALYSIS 
POSSIBLE  PARAMETERS 

Relating  model  parameters  to  physical  basin  character¬ 
istics  is  the  major  purpose  of  this  study.  The  geomorphic 
and  geophysical  structures  of  a  basin  play  an  important 
part  of  the  watershed's  response  to  rainfall.  Many  basin 
characteristics  have  been  used  in  the  past  to  relate  a 
response  of  precipitation  to  physical  structures  in  a  water¬ 
shed.  Linsley  (1982.  p.311-316)  defines  several  character¬ 
istics  of  interest  in  this  study.  They  arei  drainage  area, 
stream  density,  basin  shape,  channel  slope,  and  channel 
length . 

In  a  very  large  study  of  the  Northeastern  U.S., 

Langbein  (194?,  p.125)  using  U.S.  Geological  Survey  (USGS) 
topographic  maps  of  340  basins  withdrainage  areas  of  1.64 
to  7.797  square  miles  collected  data  on  drainage,  area, 
length  of  streams,  stream  density,  land  slope,  channel 
slope,  area-altitude  distribution,  and  area  of  water  bodies 
of  the  basins.  Characteristics  were  divided  into  geographic 
such  as  water  bodies,  direction  of  stream  flow,  latitude 
and  longitude, and  topographic  such  as  basin  area,  stream 
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length,  area-distance  distribution,  land  slope,  basin 
altitude,  and  tributary  and  principal  stream  slopes.  A 
principal  channel  was  defined  as  one  that  drains  more 
than  ten  percent  of  the  total  area  while  a  tributary 
channel  drains  less  than  ten  percent  of  the  total  area. 

Their  conclusions  were  that  many  characteristics  have 
a  direct  relationship  with  other  characteristics {  for 
example,  steep  land  slopes  imply  steep  channel  slopes 
and  stream  density  tends  to  vary  with  the  land  slope. 

The  study  showed  that  no  one  element  is  unique  in  any 
one  basin  and  that  not  all  characteristics  must  be  correla¬ 
ted  to  model  parameters  to  get  satisfactory  regression  re¬ 
lationships  . 

SOME  UNIT  HYDROGRAPH  APPLICATIONS 

Taylor  and  Schwarz  (1952,  p.235)  published  a  study 
based  on  65  rainfall  excess  periods  over  20  basins  in  the 
North  and  Middle  Atlantic  States  with  varying  drainage  areas 
of  20  to  1600  square  miles.  Their  results  were  that  the 
drainage  area,  length  of  longest  watercourse,  length  to 
center  of  area  and  the  equivalent  mainstream  slope  were 
the  most  significant  basin  characteristics. 

Brater's  (19^0,  p. 115*0  study  indicated  that  the 


unit  hydrograph  method  is  one  of  the  best  practical  devices 
for  predicting  flood  flows.  His  work  based  on  22  small 
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watersheds  in  the  Southern  Appalachians  of  varying 
cover  with  areas  of  4.24  to  1,876.7  acres  showed  that 
the  unit  hydrograph  principle  "permits  the  engineer  to 
estimate  not  merely  the  peak  discharge,  but  the  entire 
hydrograph  of  runoff.  The  peak  may  be  determined  for  any 
desired  time  interval." 

Analyzing  data  on  nine  streams  tributary  to  the 
Chemung  River  in  New  York,  Morgan  and  Hullinghorst  (1939. 
p.  1)  found  good  correlation  to  the  area  of  the  watershed, 
mean  length  of  travel  and  the  mean  height  of  the  watershed 
above  the  outflow. 

Unit  hydrographs  were  used  in  a  flood- frequency  work 
by  Kinnison  and  Colby  (1945).  On  about  48  basins  in  Mass¬ 
achusetts,  the  characteristics,  found  to  be  most  influential, 
were i  median  altitude  (s),  drainage  area  (M)  and  the  average 
distance  which  water  from  runoff  uniformly  distributed  over 
the  basin  must  travel  to  the  outlet  (L) .  Thus  M  is  a  meas¬ 
ure  of  the  volume  of  water  to  be  discharged,  s  is  a  measure 
of  the  fall  and  L  is  a  measure  of  the  distance  that  the 
water  travels  to  the  point  of  discharge. 

More  recent  works  by  Rodriguez  and  Gonzalez  (1982,  p. 
877)  and  Rodriguez  et  al.(1982,  p.  887)  have  applied  the 
geomorphoclimatic  theory  of  the  IUK.  Their  works  have 
helped  to  explain  much  of  the  noise  observed  in  relating  IUH 
parameters  to  basin  characteristics  without  considering  the 
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coupling  effects  of  the  climate  and  geomorphology  on  a 
basin. 

Hence,  it  seems  that  most  authors  use  easily  definable 
basin  characteristics  such  as  main  channel  length,  main 
channel  slope  and  drainage  area  to  obtain  their  hydrograph 
response  relationships  to  basin  charactersitics. 


BASE  FLOW 

The  release  of  water  from  underground  storage  into  the 
channel  is  called  base  flow.  Dooge  (19?3»  P*  89)  has 
presented  a  thorough  discussion  of  the  subject  as  well  as 
its  separation  from  the  total  storm  hydrograph.  Although 
several  methods  exist  to  calculate  base  flow,  the  difference 
between  the  methods  is  not  significant  in  relation  to  the 
total  flow  in  the  channel  during  a  major  storm. 


CHAPTER  THREE 
METHOD 

PURPOSE  AND  SCOPE 

The  purpose  of  this  chapter  is  to  describe  the  method 
used  in  determining  HEC-1  parameters  for  the  study  re¬ 
gion's  gaged  basins  and  in  developing  the  relationships 
between  unit  hydrograph  parameters  and  basin  characteristics. 
The  scope  of  this  chapter  is  to  present  the  HEC-1  options 
employed  and  the  multiple  regression  scheme  used  in  this 
thesis . 


WHAT  IN  THE  HECK  IS  HEC-1? 

HEC-1  is  the  shortened  name  of  HEC-1  Flood  Hydrograph 
Package,  Computer  Program  723-x6-L2010 ,  with  its  latest 
revision  taking  place  in  September  1981.  The  program 
was  developed  at  the  Hydrologic  Engineering  Center  (HEC), 
Davis,  California,  which  is  a  part  of  the  U.S.  Army 
Corps  of  Engineers'  Water  Resources  Support  Center.  The 
original  development  of  HEC-1  was  made  in  1967  by  Leo- Beard 
and  others  (HEC-1,'83),  with  the  first  version  being  pub¬ 
lished  in  October  1968.  Major  revisions  have  been  made  in 
1969»1970,  1973i  and  1981.  The  current  version  has  the 
major  additional  capabilities  of  the  dam-break  (HEC-1DB) , 
project  optimization  (HEC-1G-S),  and  the  kinematic  wave 
(HEC-1KW )  programs  included  in  the  package. 


V 
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The  HEC-1  model  is  designed  to  simulate  the  surface 
runoff  response  of  a  river  basin  to  precipitation  by 
representing  the  basin  as  an  interconnected  system  of 
hydrologic  and  hydraulic  components  within  a  portion  of 
the  basin,  commonly  called  a  subbasin.  A  component  repre¬ 
sents  a  surface  runoff  entity,  a  stream  channel,  or  a 
reservior.  Representation  of  a  component  requires  a  set 
of  parameters  which  specify  the  particular  characteristics 
of  the  component  and  mathematical  relations  which  describe 
the  physical  processes.  The  result  of  the  modelling  process 
is  the  computation  of  streamflow  hydrographs  at  desired 
locations  in  the  river  basin. 

As  previously  stated  in  the  Introduction  under  Scope 
and  Goals,  this  thesis  will  use  HEC-1  to  model  gaged  basins 
incrder  to  determine  values  for  parameters  of  a  surface 
runoff  entity  component  which  will  be  used  as  a  predicting 
tool  for  ungaged  basins  in  the  region  of  interest 

ASSUMPTIONS  AND  LIMITATIONS 

As  pointed  out  in  the  HEC-1  Users  Manual,  HEC-1  makes 
a  number  of  assumptions  and  has  some  limitations.  The  major 
assumptions  are  that  hydrologic  processes  such  as  precip¬ 
itation  and  interception/infiltration,  may  be  represented 
by  model  parameters  which  reflect  average  conditions  within 
a  subbasin j  model  parameters  represent  temporal  ss  well  as 
spatial  averages.  Some  of  its  limitations  are  that  simu- 


lation  is  limited  to  a  single  storm  event  (because  HEC-1 
does  not  have  the  ability  to  account  for  soil  moisture 
recovery  during  periods  of  no  precipitation)  and  that 
results  are  given  in  discharge  rates  rather  than  in  stage 
heights. 

HYDROLOGIC  ANALYSIS  PROCEDURES  OF  UNGAGED  WATERSHEDS  USING 

HEC-1 

Guidelines  and  methods  for  modelling  ungaged  basins 
are  presented  in  a  HEC  report  entitled  "Hydrologic  Analysi 
of  Ungaged  Watersheds  Using  HEC-1 ",  Training  Document  No. 
15.  April  1982.  This  report  formed  the  methodology  used 
in  this  thesis.  The  report  gave  detailed  techniques 
for  estimating  and  calibration  of  HEC-1  model  parameters 
including  discussions  of  the  runoff  transformation  and 
loss  rate  parameters. 

Below  is  outlined  the  basic  method  used  in  the 
regional  analysis  of  watershed  characteristics  of  this 
thesis  (which  is  described  in  greater  detail  in  the 
Report) » 

1.  Collect  precipitation  and  runoff  information 
for  a  range  of  major  flood  events  (usually  5-iO  for 
each  gaged  basin)  in  the  region.  Identify  the 
watershed  boundaries  on  a  topographic  map  where 
the  gaging  stations  are  taken  as  the  basin  outlets. 

2.  Run  a  HEC-1  rainfall-runoff  analysis  to  optimize 
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the  unit  hydrograph  and  loss  rate  parameters  for 
the  gaged  drainage  areas. 

3.  Correlate  the  unit  hydrograph  and  loss  rate 
parameters  with  basin  characteristics  and  develop 
generalized  relationships  for  these  parameters. 

4.  Use  the  generalized  relationships  to  compute 
parameters  for  the  ungaged  basins  by  means 

of  measurable  basin  characteristics. 

5.  Perform  a  watershed  simulation  using  HEC-1 
for  several  historical  storm  events  for 
verification  of  reconstituted  runoff  peaks  and 
volumes  at  gaged  locations,  if  unsatisfactory 
results  are  obtained,  the  analysis  is 
repeated  after  adjusting  the  parameters. 


RAINFALL-RUNOFF  SIMULATION 


The  rainfall-runoff  simulation  involves  five  processes. 


They  are  precipitation,  intercept ion/ infiltration,  trans¬ 
formation  of  precipitation  excess  to  basin  outflow,  base- 
flow,  and  flood  routing  (not  used  in  this  study).  Each  of 
the  first  four  processes  will  be  discussed  below. 


PRECIPITATION 

The  precipitation  distribution  chosen  for  this 
study  among  the  possibilities  in  HEC-1  was  the  method 
of  weighted  precipitation  gages.  This  method  determines 
a  basin  average  rainfall  amount  for  a  historical  storm 
and  distributes  it  temporally  over  the  entire  basin. 

This  process  produces  a  precipitation  hyetograph. 

Each  gaged  basin  in  the  study  region  was  divided 
into  Thiessen  polygons  (see  Linsley  (1982,  p.  71)  for 
the  Thiessen  polygon  method).  With  one  rain  gage  per 
polygon,  the  area  of  each  polygon  covering  the  basin 
could  be  determined  using  a  planimeter.  The  area  of  a 
polygon  divided  by  the  total  basin  drainage  area  yielded 
the  percent  (expressed  as  a  decimal)  of  the  relative 
weight  for  each  rain  gage.  Thiessen  polygons  were  drawn 
for  the  recording  (hourly  measurements)  and  non-recording 
(daily  measurements)  rain  gages.  To  obtain  the  total 
storm  precipitation  the  relative  weights  were  placed 
on  the  PT-FW  input  cards  (input  card  data/format  will 
be  discussed  later  in  Chapter  Four)  for  the  recording 
and  non-recording  gages  which  actually  contributed  rain 
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to  the  basin.  However,  in  order  to  develop  the  temporal 
pattern  on  the  basin, the  relative  weights  placed  on  the 

PR-PW  input  cards  only  used  recording  rain  gage  polygons. 

The  total  storm  precipitation  for  a  basin  was  computed 
as  the  weighted  average  of  measurements  from  several  gages 
using  the  formula*  h 


PRCPN(J)  *  WTN(J) 


PRCPA 


J=i 

n 

^  WTN(J) 


where  PRCPA  is  the  basin  average  total  precipitation, 
PRCPN(J)  the  total  precipitation  for  gage  J,  WTN(J)  the 
relative  weight  for  gage  J,  and  n  the  number  of  gages. 

The  temporal  pattern  for  distribution  of  the  storm- 
total  precipitation  is  computed  as  a  weighted  average  of 
temporal  distributions  from  recording  stations  using 
the  formula «  p 

2  PRCPR ( I , J )  *  WTR(J) 

prcp(i)  =  - sii! - 

\  WTR(J) 

J=> 

where  PRCP(I)  is  the  basin-average  precipitation  for  the 
Ith  time  interval,  PRCPR(I,J)  the  recording  station  precip¬ 
itation  for  the  gage  J. 

The  basin-average  precipitation  hyetograph  is  then 
computed  using  the  temporal  pattern,  PROP,  to  distribute 
the  total  rainfall  amount,  PRCPA,  see  Figure  3.1. 


to  tj  tt  t,  t4  t, 

TIME  (hr) 

Figure  3.1.  Precipitation  Hyetograph  (HEC-1,  1981). 

INTERCEPTION/INFILTRATION 

Definitions.  Interception  and  depression  storage 

are  intended  to  represent  the  surface  storage  of  water 

by  trees  or  grass,  local  depressions  in  the  ground  surface, 

in  cracks  and  crevices  in  parking  lots  or  roofs,  or  in  a 

surface  area  where  water  is  not  free  to  move  as  overland 

flow.  Infiltration  represents  the  movement  of  water  to 

areas  beneath  the  land  surface.  Both  of  these  are  lumped 

together  in  HEC-1  and  considered  as  precipitation  losses. 

Two  important  factors  should  be  noted  about  this 

process  1  one,  precipitation  which  does  not  contribute  to  'the 
runoff  process  is  considered  to  be  lost  from  the  system 

and  two,  there  is  no  provision  in  the  model  for  soil  mois¬ 
ture  or  surface  storage  recovery.  The  second  factor  has 
already  been  mentioned  under  the  Assumptions  and  Limitations 
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portion  of  this  chapter. 

As  in  the  precipitation  process,  the  interception/ 
infiltration  process  uniformly  distributes  the  precipi¬ 
tation  loss  over  the  entire  basin.  This  process  yields 
a  basin-average  for  rainfall  losses.  There  is  a  feature 
in  the  model  in  which  a  percentage  of  the  basin  may  be 
labeled  as  impervious  and  no  losses  are  calculated  for  this 
portion;  however,  this  feature  was  not  used  in  this  study. 

Four  different  methods  may  be  used  in  calculating 
the  precipitation  loss.  By  using  any  one  of  these  methods, 
an  average  precipitation  loss  is  determined  for  a  computa¬ 
tion  interval  and  subtracted  from  the  rainfall  hyetograph. 
The  resulting  precipitation  excess  is  used  to  compute 

an  outflow  hydrograph  for  a  subbasin  as  is  shown  in 


Figure  3.2. 
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Input  Hyetograph  (HEC-1,  1981). 


The  four  methods  available  to  calculate  rainfall 
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losses  are  initial  and  uniform  loss  rate,  exponential 
loss  rate,  Soil  Conservation  Service  (SCS)  curve  number, 
and  the  Holtan  loss  rate.  This  thesis  only  considered  the 
initial  and  uniform  loss  rate.  It  will  be  the  only  method 
discussed  here. 

The  initial  and  uniform  loss  rate  method  has  two  par¬ 
ameters;  an  initial  loss,  STRTL  (units  of  depth),  and  a 
constant  loss  rate,  CNSTL  (units  of  depth/hour).  All  rain¬ 
fall  is  lost  until  the  volume  of  initial  loss  is  satisfied. 
After  the  initial  loss  is  satisfied,  rainfall  is  lost  at 
the  constant  rate,  CNSTL.  These  parameter  values  are  entered 
as  input  on  the  LU  card. 


TRANSFORMATION  OF  PRECIPITATION  EXCESS  TO  BASIN  OUTFLOW 

The  transformation  of  rainfall  excess  into  basin  out¬ 
flow  for  this  thesis  is  accomplished  through  the  unit 
hydrograph  technique,  more  specifically,  the  technique 
that  was  proposed  by  Clark  (19^5)  which  has  already  been 
discussed  in  Chapter  Two . 

The  basic  methodology  of  any  unit  hydrograph  component 
inHEC-1  is  the  same.  The  rainfall  excess  hyetograph  is 
transformed  to  a  subbasin  outflow  by  utilizing  the  general 
equation;  N  i 

Y  Y 

Q(i)  »  £i  £1  OU)  *  *<i-d+l) 

where  Q(i)  is  the  subbasin  outflow  at  the  end  of  computation- 


al  interval  i,  U ( j )  the  jth  ordinate  of  the  unit  hydrograph, 
X(i)  the  average  rainfall  excess  for  the  computational 

interval  i,  and  N  the  number  of  rainfall  ordinates. 

Application  of  the  unit  hydrograph  technique  implies 
two  assumptions*  one,  the  unit  hydrograph  is  characteristic 
of  the  subbasin  and  is  not  storm  dependent  and  two,  runoff 
due  to  excess  from  different  periods  of  rainfall  excess  may 
be  linearly  superposed. 

The  Clark  unit  hydrograph  requires  three  parameters 
to  calculate  a  unit  hydrograph j  the  time  of  concentration 
for  the  subbasin,  TC  (units  of  time),  a  storage  coefficient, 

R  (units  of  time),  and  a  time-area  curve. 

The  time  of  concentration,  TC,  may  be  defined  conceptu¬ 
ally  as  the  amount  of  time  required  for  a  particle  of 
water  deposited  at  the  furtherest  point  upstream  in  the 
watershed  to  travel  through  the  watershed  until  it  passes 
the  subbasin  outlet.  TC  is  estimated  by  the  lag  time  between 
the  end  of  the  runoff  producing  rainfall  to  the  inflection 
point  on  the  recession  limb  of  the  surface  runoff  hydrograph. 

The  storage  coefficient,  R,  is  used  to  account  for 

the  effect  of  subbasin  storage  on  the  hydrograph.  This 
parameter  is  estimated  by  dividing  the  flow  at  the  recession 
inflection  point  of  the  surface  runoff  hydrograph  by  the 
rate  of  change  of  discharge  (slope)  at  this  same  time. 

See  Figure  3.3. 
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The  time-area  curve  defines  the  cumulative  area  of  the 
watershed  contributing  runoff  to  the  subbasin  outlet  as  a 
function  of  time  (expressed  as  a  portion  of  TC). 

No  developed  time-area  curves  for  any  of  the  gaged 
subbasins  in  the  study  region  could  be  found,  therefore, 
use  of  the  HEC-1  synthetic  dimensionless  time-area  curve 
was  made.  The  equations  of  this  curve  are t 

AI  =  1.414  T1*5  0  *  T*.5 

1  -  AI  =  .1.414  (1-T)1’5  .5  -  T  -  1 


where  AI  is  the  cumulative  area  as  a  fraction  of  the  total 
subbasin  area  and  T  the  fraction  of  the  time  of  concentra¬ 
tion.  See  figure  3.4. 
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The  ordinates  of  the  time-area  curve  are  converted 
to  volume  of  runoff  per  second  for  unit  excesr.  and  inter¬ 
polated  to  the  given  time  interval.  This  results  in  the 
development  of  a  translation  hydrograph  which  is  then 
routed  through  a  linear  reservoir  to  simulate  the  storage 
effects  of  the  subbasin.  The  routing  produces  a  unit  hydro 
graph  for  instantaneous  excess  which  is  averaged  to  yield 
the  hydrograph  for  unit  excess  occurring  in  the  given 
time  interval  (one  hour  in  all  cases  of  this  thesis). 

The  linear  reservoir  routing  is  accomplished  using 
the  equation 

Q(2)  »  CA  *  I  ♦  CB  *  Q(l) 
where  the  routing  coefficients  are 
CA  =  ot/(R  +  .5at) 

CB  =  1  -  CA 


QUNGR  =  .5(  Q(1 )  ♦  Q(2)  ) 

and  where  Q(2)  is  the  instantaneous  flow  at  the  end  of 
the  period,  Q(l)  the  instantaneous  flow  at  the  beginning 
of  the  period,  I  the  ordinate  of  the  translation  hydrograph, 
At  the  computational  time  interval  (one  hour  in  all  cases 
in  this  study),  R  as  previously  defined,  and  QUNGR  the  unit 
hydrograph  at  the  end  of  the  computation  interval. 

The  parameters  TC  and  R  are  entered  as  input  on  the 
UC  card. 

BASEFLOW 

Unlike  direct  surface  runoff  which  is  calculated 
as  the  total  precipitation  minus  the  losses,  baseflow 
is  defined  as  the  release  of  water  from  subsurface  storage. 
To  include  baseflow  effects  the  HEC-1  model  must  be  supplied 
with  three  input  parameters, STRTQ,  QRCSN,  and  RTIOR. 

STRTQ  is  the  initial  discharge  (starting  flow)  in 
the  stream  outlet  at  the  beginning  of  the  storm.  It  is 
the  flow  at  the  outlet  at  the  beginning  of  the  rainfall 
hyetograph.  It  is  affected  by  the  long  term  contribution 
of  groundwater  releases  in  the  absence  of  precipitation  and 
is  a  function  of  antecedent  conditions. 

QRCSN  is  the  flow  on  the  receding  limb  of  the 
computed  hydrograph  at  which  an  exponential  recession 
begins. 

RTIOR  is  a  user  specified  exponential  decay  rate 


which  is  assumed  to  be  characteristic  of  the  subbasin. 

It  is  equal  to  the  ratio  of  a  recession  limb  flow  occurring 
one  hour  later.  See  Figure  3.5. 


Figure  3.5*  Base  Flow  Parameters  (HEC-1,  1981). 


Baseflow  or  recession  flow  is  computed  with  the 
equation  1 

Q  =  Qo  <RTIOR)n 

where  Q  is  STRTQ  or  QRCSN  and  n  the  number  of  time  intervals 
since  recession  was  initiated. 

The  values  of  QRCSN  and  RTIOR  are  obtained  by  plotting 
the  observed  flows  (starting  at  the  peak)  versus  time  on 
semi-log  paper.  The  point  at  which  the  flows  begin  to 
become  a  straight  line  defines  QRCSN.  See  Figure  3.6. 


Log  Q 


RTIOR  =  Qa/Qb 


RTIOR  is  usually  obtained  by  taking  the  average 
ratio  of  Qa  and  Qt,  (flows  one  hour  apart)  of  several, 
usually  about  ten,  along  the  straight  line  portions  of  Figure 
3.6  In  this  study  a  subbasin  average  is  determined  by 
taking  the  average  RTIOR  for  each  storm  on  the  basin  and 
then  averaging  these  values.  The  subbasin  average  was 
then  held  constant  throughout  all  calculations  for  that 
subbasin  while  using  HEC-1. 

The  computed  flood  hydrograph  is  adjusted  to  include 
the  baseflow.  At  the  beginning  of  a  storm  the  baseflow 
is  computed  starting  at  STRTQ  and  decays  according  to  the 
formula  above  until  the  computed  flood  flow  on  the  falling 
limb  equals  the  value  of  QRCSN.  From  then  on,  the  flood 
flow  is  calculated  using  the  base  flow  equation  as  the 
computed  flood  flow  unless  the  computed  flow  rises  above 
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the  recession  flow.  This  would  occur  in  the  case  of  a 
double-peaked  hydrograph  such  as  in  Figure  3.7. 


CALIBRATION  OF  MODEL  PARAMETERS 
As  with  all  numerical  models,  the  HEC-1  rainfall- 
runoff  model  for  single  event  simulation  requires  data 
for  calibration  in  order  to  be  a  predictive  tool.  The 
parameters  of  the  numerical  model  are  determined  by  using 
observed  precipitation  and  runoff  data  to  solve  the 
inverse  problem?  in  other  words,  given  the  system  input 
(precipitation)  and  system  output  (runoff),  the  inverse 
problem  consists  of  defining  the  characteristics  of  a 
system,  that  produces  the  transformation  from  input  to 
output.  HEC-1  has  two  approaches  to  transform  rainfall 
values  to  runoff.  They  are  the  unit  hydrograph  and 
kinematic  wave  approaches .  The  unit  hydrograph  is  the 
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is  the  most  commonly  used  and  was  the  approach  selected 
for  this  study.  The  unit  hydrograph  approach  assumes 

that  a  single  unit  hydrograph  is  appropriate  for  all 
magnitudes  of  rainfall  excess.  HEC-1  has  three  unit 
hydrograph  methods  available  in  its  package,  those  of 
Clark,  Snyder,  and  the  Soil  Conservation  Service  (SCS). 

The  U.S.  Army  Corps  of  Engineers  uses  the  Clark  unit 
hydrograph  technique  most  frequently  and  for  this  reason 
was  the  technique  selected  for  this  thesis. 

HEC-1  has  the  capability  of  automatically  producing 
its  "best"  estimate  of  selected  model  parameters  through 
its  optimization  subroutine.  This  subroutine  selects  those 
values  of  the  parameters,  which  yield  the  "best"  repro¬ 
duction  of  some  measured  runoff  event  with  the  available 
measured  precipitation  data  and  the  selected  modelling 
approach.  This  automatic  calibration  approach  is  accom¬ 
plished  through  an  objective  function,  defined  asj 


STDER  = 


(  (QOBSi  “  QCOMPi)2 


WTi  ) 


where  STDER  is  the  error  index,  QOBSi  "the  observed  runoff 
hydrograph  ordinate  for  period  i,  QCOMPi  the  computed 
runoff  hydrograph  ordinate  for  period  i,  computed  by 
HEC-1  with  current  parameter  estimates.  N  is  the  total 
number  of  hydrograph  ordinates  and  WTi  a  weight  for  the 
hydrograph  ordinate  defined  asi 


WTi  =  (QOBSi  ♦  QAVE)/(2  *  QAVE) 
where  QAVE  is  the  average  computed  discharge. 

The  5TDER  calculation  for  optimization  may  be  viewed 
graphically  as  in  Figure  3.8, 


Figure  3.8.  Optimization  Calculation  (HEC-1,  1981) 

in  which  the  limits  of  optimization  may  be  user  specified 
on  the  OU  input  card  in  order  to  get  a  better  fit  on  the 
portion  of  the  curve  which  is  of  interest. 

The  WTi  'tenn  biases  the  objective  function  by  em¬ 
phasizing  accurate  reproduction  of  peak  flow  rather  than 
low  flow.  The  objective  of  this  function  is  to  minimize 
the  value  of  STDER  in  order  to  get  the  "best'’  fit.  The 
word  "best"  has  been  placed  in  quotation  marks  to  indicate 
that  the  values  determined  through  this  optimization 
technique  do  not  guarantee  that  the  values  are  global 
solutions;  indeed  an  example  of  this  is  shown  in  Chapter 
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Five.  A  better  solution  may  be  found  by  varying  the 
starting  values  of  the  parameters  from  their  default  values 
selected  by  HEC-1  in  its  first  run.  These  default  values 
are  listed  below t 

Parameter  Name 
TC  ♦  R 
V(TC  +  R) 

STRTL 
CNSTL 


.2, 


Initial  Value 
(Drainage  Area  (mi*)) 
0.5 
1.00 
0.1 


1/2 


The  user  also  has  the  option  of  designating  initial 
values  for  model  optimization  as  well  as  designating 
which  parameters  are  to  be  held  constant  during  the  opti¬ 
mization.  Assignment  of  the  initial  values  is  discussed  in 
the  HEC-1  Programmer’s  Manual,  October  1981.  Constraints 
are  placed  on  parameter  values.  These  values  are  bounded 
by  a  range  of  feasible  values  because  of  the  physical 
limitations  on  the  values  that  the  various  unit  hydrograph 
and  loss-rate  parameters  may  have.  These  constraints  are» 

TC  +  R  -  I.03  At/(1-R/(TC+R)) 

R  *  .52 At 

where  t  is  the  computation  interval  which  is  one  hour  in 
all  cases  in  this  study. 

STRTL  -  0 
CNSTL  -  0 


A  complete  discussion  of  the  optimization  technique 
as  well  as  the  algorithm  is  presented  in  the  HEC-1 
User's  Manual.  The  technique  is  a  univariate  search 
technique  in  which  each  parameter  value  to  be  optimized 
is  varied  one  at  a  time  with  all  other  parameters  held 
constant.  The  "best"  value  of  each  parameter  is  estimated 
hy  Newton's  method.  After  each  parameter  is  estimated 
in  turn  through  four  complete  cycles,  the  parameter 
which  most  improved  the  objective  function  in  its  last 
change  is  adjusted  again.  This  adjustment  continues  until 
no  one  parameter  yields  a  reduction  in  the  objective 
function  of  more  than  one  percent.  After  one  more  complete 
search  of  all  parameters,  a  final  adjustment  to  the  computed 
hydrograph  volume  to  within  one  percent  of  the  observed 
hydrograph  volume  is  made.  The  final  objective  function 
value  is  determined  from  this  final  adjustment .  An  example 
of  the  optimization  output  is  shown  in  Figure  3.9. 

Listed  first  are  the  initial  values  of  the  parameters. 
The  asterisked  (*)  values  denote  which  variable  was 
changed  and  its  "optimum"  value  along  with  the  objective 

function  value  with  the  other  parameters. 

After  the  objective  function  routine  is  printed,  the 
optimization  results  are  blocked  out  followed  by  a 
statistical  summary.  If  the  user  specifies  optimization 
limits  on  the  OU  card,  two  summaries  are  printed.  The 
first  contains  statistics  based  on  the  optimization  region 
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while  the  second  summary  is  based  upon  the  entire  flood 
hydrograph . 


The  definition  of  the  terms  used  in  the  summary  are 
self-explanatory  except  for  the  following  four  terms; 
Standard  Error- the  root  mean  squared  sum  of  the 
difference  between  observed  and  computed  hydrographs. 
Objective  Function — the  weighted  root  mean  squared 
sum  of  the  difference  between  observed  and  computed 
hydrographs . 

Average  Absolute  Error — the  average  of  the  absolute 
value  of  the  differences  between  observed  and 
computed  hydrographs. 

Average  Percent  Absolute  Error — the  average  of 
absolute  value  of  percent  difference  between  observed 
and  computed  hydrograph  ordinates. 

APPLICATION  OF  THE  CALIBRATION  CAPABILITY 
The  following  steps  represent  the  strategy  recom¬ 
mended  to  determine  the  model  parameters  for  HEC-1  in 
this  study; 

1.  Por  each  storm!  the  base  flow  and  recession 
parameters  (STRTQ,  QRCSN ,  and  RTIOR)  are  determined 
graphically.  These  parameters  can  not  be  estimated 
automatically  and  must  be  entered  as  input.  These 

parameters  are  placed  on  the  BF  input  card  for  each 
storm. 
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2.  For  each  storm  and  each  gage,  the  optimal  estimates 
of  all  unknown  unit  hydrograph  and  loss  rate  parameters 
(TC,  R,  STRTL,  and  CNSTL)  are  determined  using  the 
automatic  calibration  feature. 

3.  A  value  of  the  constant  loss  rate,  CNSTL,  is  esti¬ 
mated.  A  regional  value  of  CNSTL  based  on  all  storms 
at  all  the  gages  is  selected. 

4.  With  CNSTL  fixed  at  the  selected  value,  the 
automatic  calibration  feature  is  again  employed  for 
all  the  storms  at  all  the  gages.  A  regional  relation¬ 
ship  for  the  initial  loss,  STRTL,  is  selected.  (See 
Figure  5-1  ) 

3.  With  both  loss  rate  parameters  fixed  the  automatic 
calibration  of  the  model  is  repeated. 

6.  Values  of  TC  +  R  are  selected  for  each  basin. 

An  average  value  of  R/(TC  +  R)  is  selected  for  the 
region . 

7-  After  all  parameters  are  selected,  the  values  are 
verified  by  simulating  the  response  of  the  gaged 
basins  to  other  storms  not  used  in  the  calibration 
process . 

REGRESSION  ANALYSIS  APPROACH  TO  PARAMETER  ESTIMATION 

(HEC-1,  1982) 

Multiple  linear  regression  techniques  were  used  to 
correlate  the  unit  hydrograph  parameters  with  several 
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basin  characteristics.  Nonlinear  relationships  were  inves¬ 
tigated  by  transforming  values  logarithmically.  Multiple 
regression  analysis  requires  a  computer  program.  Herein, 
the  computer  package  MINITAB  was  used  in  the  analysis. 

The  parameters  TC  and  R  were  considered  dependent 
variables  and  the  basin  characteristics  were  the  inde¬ 
pendent  variables.  TC,  R,  (TC  ♦  R),  and  R/(TC  +  R)  were 
analyzed,  first,  by  considering  several  basin  character¬ 
istics  and,  later,  by  reducing  the  number  of  independent 
variables  until  the  "best"  relationships  (equations) 
were  found. 

The  statistics  describing  the  "goodness-of-fit" 
of  the  regression  equation  to  the  data  are  used  in  evalu¬ 
ating  the  analysis .  These  statistics  are  the  coefficients 
of  determination  (both  the  adjusted  and  unadjusted), 

the  partial  determination  coefficient,  and  the  standard 
error  of  estimate.  Their  definitions  are  given  belowi 

1.  The  adjusted  and  unadjusted  multiple-determination 

2 

coefficients  (R  )  are  a  measure  of  the  percent  of 
variance  in  the  dependent  variable  explained  by 
the  independent  variable.  The  magnitude  of  these 
coefficients  varies  from  0  to  1.  The  closer  to 
the  value  of  1,  the  greater  the  reliability  is 
of  the  estimate . 

2,  The  partial-determination  coefficient  (  r2)  is 

a  measure  of  the  importance  of  an  independent  variable 
by  determining  the  reduction  in  variance  in  the 
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dependent  variable  when  the  variable  is  included 
with  the  other  independent  variables. 

3*  The  standard  error  of  estimate  (3e)  is  the  standard 
deviation  of  the  difference  between  the  observed 
dependent  values  and  the  values  computed  from  the 

regression  equation  in  the  units  of  the  dependent 

variable;  therefore,  it  must  be  compared  with  the 

mean  and  the  standard  deviation  of  that  variable. 

For  a  more  detailed  description  of  multiple-regression 

analysis  the  reader  is  referred  to  a  statistics  book,  such 

as  Benjamin  (1970,  p.4l9).  The  general  rule  is  to  use 
o  2 

the  values  of  R  ,  R  adjusted,  and  Se  computed  for  each 

-  ’  .j- 

regression  equation  as  a  guide  and  select  the  equation 

with  the  fewest  independent  variables  and  the  best  values 
of  R2  and  Se* 


CHAPTER  FOUR 
DATA 

PURPOSE  AND  SCOPE 

The  purpose  of  this  chapter  is  to  present  the  data 
used  in  this  study.  The  scope  of  this  chapter  is  to  dis¬ 
play  the  basic  data  obtained  from  various  sources  after 
some  refinement  into  a  usable  form  for  this  analysis. 

BASIN  CHARACTERISTICS 

Nine  watersheds  in  the  Susquehanna  and  Chemung  River 
basins  were  chosen  for  the  present  investigation.  The 
watersheds  were  unregulated  and  possessed  a  U.S.  Geological 
Survey  (U3GS )  stream  gage.  The  locations  of  these  gages 
were  taken  as  the  outlets  from  the  watersheds .  The  nine 
basins  are  henceforth  referred  to  by  the  main  channel 
watercourse.  Their  names  arei  Butternut  Creek,  Canasawacta 
Creek,  Charlotte  Creek,  Corey  Creek,  Elk  Run  Creek,  Five 
Mile  Creek,  Newtown  Creek,  Otselic  River,  and  Owego 
Creek.  All  basins  are  in  New  York  state  except  Corey 
Creek  and  Elk  Run  Creek,  which  are  in  Pennsylvania.  Figure 
4.1  depicts  the  basins  used  and  their  location  with  re¬ 
spect  to  one  another. 

The  basin  characteristics  were  measured  from  USGS 
Topographic  Quadrangle  maps  (  1»2400  ).  Table  4.1  contains 
a.  tabular  synopsis  of  the  basin  characteristics  chosen 
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for  the  present  analysis  on  the  basis  of  the  review  in 
Chapter  2.  A  definition  of  each  term  follows »  (after  Linsley 
(1955.  p-  313)) 

DRAINAGE  AREA  (DA) — The  area  in  square  miles  of  the  water¬ 
shed  above  the  stream  gage.  Any  excess  rainfall  falling 
in  the  watershed  would  at  some  time  pass  through  the 
gaging  station.  Rainfall  falling  outside  of  this  area 
would  drain  into  another  watershed.  These  data  were  ob¬ 
tained  from  USGS  Water  Supply  Papers. 

SLOPE  (S) — This  is  the  average  main  channel  slope  in 
vertical  feet  per  horizontal  channel  mile.  Slopes  were 
obtained  by  plotting  the  elevation  of  the  map  contours 
(in  feet)  versus  the  distance  of  the  map  contours  from 
the  gaging  stations  (in  miles)  and  taking  the  straight 
line  slope  fitted  by  the  least  squares  method.  The  steeper 
the  slope  the  faster  the  basin  will  drain. 

SLOPE  2(S2) — This  is  similar  to  slope  (S)  except  the  main 
channel  is  divided  into  two  segments.  Two  straight  lines 
are  drawn  representing  a  closer  fit  to  the  actual  main 
charnel  profile.  Slope  2  is  a  weighted  average  of  these 
two  segments. 

S2  =( (slope  of  segment  1 ) (horizontal  length  of 
segment  l)  ♦  (slope  of  segment  2) (horizontal 
length  of  segment  2))/(total  main  channel  length) 


SLOPE  WT  (SWT) — This  slope  is  a  weighted  slope  intended  to 
represent  a  better  picture  of  the  basin's  slope  character¬ 
istics  by  including  the  slope  contribution  of  some  of  the 
tributary  streams.  Its  feet  is  also  listed  in  feet  per 
mile . 

DRAINAGE  DENSITY  (DD) — This  is  a  measure  of  the  amount 
of  stream  channel  miles  per  basin  square  mile.  It  is 
obtained  by  measuring  the  cumulative  length  of  every 
stream  channel  shown  on  the  USGS  topographic  map.  This 
cumulative  length  is  divided  by  the  drainage  area.  The 
higher  the  drainage  density,  the  smaller  distance  any 
particle  of  excess  rainfall  must  travel  before  entering 
a  stream  channel.  This  implies  that  the  larger  the  DD 
value,  the  faster  the  basin  will  drain. 

LENGTH  OF  THE  MAIN  CHANNEL  (L)— This  is  a  measure  of  the 
longest  watercourse  measured  from  the  gaging  station  to 
the  furtherest  point  of  stream  channel  origination  upstream. 


L  is  measured  in  miles. 

LENGTH  TO  THE  BASIN  CENTER  OF  MASS  (LCM)~ The  center  of 
mass  is  found  by  hanging  a  paper  cut-out  diagram  of  the 
drainage  basin  and  marking  off  the  vertical  lines  obtained 
from  three  different  suspension  points.  The  value  for  LCM 
is  obtained  by  measuring  the  distance  from  the  gaging 
station  to  the  point  on  the  main  channel  opposite  the 
center  of  mass.  LCM  is  measured  in  miles. 

SHAPE  FACTOR  ONE  (SF1) — This  is  a  measure  of  the  basin's 
shape.  It  is  equal  to  the  length  of  the  main  channel  squared 
divided  by  the  drainage  area,  i.e.  L  /DA. 

SHAPE  FACTOR  TWO  (SF2)— It  is  similar  to  SF1.  It  is  the 
ratio  between  the  diameter  of  a  circle,  the  size  of  the 
drainage  area  (DA),  and  the  length  of^he  main  channel  (L), 
i.e.  SF2  =  DIA/L  where  DIA  =  (4*DA/rr)  •  This  term  was 

proposed  by  Cruff  and  Rantz  (1965). 

ALTITUDE  (ALT) — This  is  a  measure  of  the  vertical  drop 
between  the  elevation  at  the  main  channel  headwaters  and 
the  gaging  station  in  feet. 

STOBM/FLOOD  DATA  SELECTION 

An  effort  was  made  to  obtain  5-10  of  the  major  flood 
events  occurring  in  each  basin  since  the  record  flooding 
which  was  produced  by  Hurricane  Agnes  (June  1972).  There 
were  several  seasons  for  selecting  events  since  the  summer 
of  1972.  First,  it  is  desired  to  minimize  the  influences 
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of  basin  changes  which  either  naturally  occur  or  are 
caused  by  man  over  time.  It  was  felt  that  a  time  period 
of  ten  to  eleven  years  would  be  a  reasonable  length  of 
time  to  minimize  basin  change  effects.  Second,  results 
of  studies  prior  to, or  as  a  result  of.  Hurricane  Agnes 
may  be  compared  with  the  results  of  this  work.  The  storm 
events  for  each  basin  are  presented  in  Table  4.2.  It 
should  be  noted  that  not  all  storm  dates  were  taken  after 
June  1972,  because  not  enough  major  storms  had  occurred 
to  provide  a  sufficiently  wide  data  base  to  conduct  the 
analysis.  In  addition,  only  events  were  selected  which  were 
free  from  any  indication  of  snowmelt,  ice  jam,  or  difficul¬ 
ties  involving  the  stream  gaging  stations.  These  problems 
are  beyond  the  scope  of  this  work. 

PRECIPITATION  DATA 

Precipitation  data  for  the  storm  events  in  each  basin 
were  obtained  from  documents  published  by  the  National 
Weather  Service.  Basins  were  traced  on  paper  from  small 
scale  maps.  Locations  of  the  rain  gages  in  the  vicinity 
of  each  basin  were  plotted  on  the  tracing  paper.  The  Thiessen 
Polygon  method  was  used  to  divide  up  the  basin.  Planimeter 
measurements  were  made  for  all  recording  and  nonrecording 
gages  influencing  the  basin.  The  fraction  of  the  basin 
contained  in  each  polygon  was  recorded  along  with  the 
rainfall  amounts  obtained  from  the  Climatological  Data 


Butternut  Creek  (NY) 
10-18-75* 
10-09-76 
10-21-76 
9-26-77 


Five  Mile  Creek 
5-07-75 
9-26-75 
9-14-77 
9-25-77  * 


10-17-77 


Newtown  Creek  (NY) 


Canasawacta  Creek  (NY) 

9-26-75 

6-26-68 

10-09-76  * 

11-18-68 

10-21-76 

6-24-69 

9-25-77 

7-31-71  * 

10-17-77 

6- 22/23-72 

7- 03-74 

11-04-77 

11-08-7? 

11-11-77 

Charlotte  Creek  (NY) 

10-06-79 

6-29-73 

11-26-79 

12-21-73 

10-28-81  * 

Corey  Creek  (PA) 

Otselic  River  (NY) 

11-14-72 

7-03-74 

9-26-75 

9-26-75 

10-09-76  * 

10-09-76  * 

5-14-78 

10-17-77 

10- 05-79 

11- 26-79 

4-05-78 

6-05-82 

Owego  Creek  (NY) 
11-08-72 

Elk  Run  Creek  (NY) 

9-26-75 

10-22-70 

10-09-76  * 

11-14-72 

9-25-77 

9-26-75 

10-09-76 

10-20-76 

5-14-78 

8-06-78  * 

10-17-77 

Total  Storms 
52 

Table  4.2 

Storm  Data  (*) — Verification  storms  not  used  in  model 

calibration . 
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and  Hourly  Precipitation  Data  (recording  stations),  docu¬ 
ments  of  the  National  Weather  Service .  Appendix  A  contains 
the  watershed  rain  gage  locations  and  the  Thiessen  polygons. 


STREAMFLOW  DATA 

The  streamflow  data  for  each  flood  event  at  a  partic¬ 
ular  gage  were  obtained  through  the  Regional  Office  of  the 

U.S.G.S.  located  in  Ithaca,  New  York  and  the  State  office 
of  the  U.S.G.S.  located  in  Harrisburg,  Pennsylvania  for 

the  Pennsylvania  stations.  Stage  data  were  obtained  as 
height  of  the  stream  at  particular  times.  For  this  study 
therefore,  the  stage  data  (FT)  had  to  be  converted  to 
flow  rates  (FT^/sec)  by  the  use  of  the  stream's  rating 
curve.  Flow  rates  were  tabulated  and  retained  for  use  as 
input  for  the  HEC-1  calibration  sequence. 

The  recession  parameters,  RTIOR  and  QRCSN  were  calcu¬ 
lated  as  described  in  Chapter  Three  and  retained  as  input 
for  the  computer  model. 

INPUT  FORMATTING  AND  DATA  STRUCTURE  FOR  HEC-1 

Input  for  the  HEC-1  model  is  accomplished  by  creating 
a  card  deck  organized  in  a  certain  structure  and  sequence 
for  easy  computer  reading.  Each  card  of  eighty  columns 
is  set  up  in  the  following  manner;  columns  1  and  2  are 
used  to  identify  the  card  type  and  purpose.  Following 
columns  1  and  2  are  ten  data  fields  of  eight  columns  each, 


except  the  first  field  which  is  composed  of  only  six  columns. 
More  detailed  discussion  of  data  organization  may  be  found 
in  HEC-1  manuals.  A  brief  line-by-line  explanation  of  a 
typical  card  deck  as  used  in  this  study  may  be  found  in 
Appendix  B.  as  well  as  the  input  card  data  information 
for  all  the  storms  used  in  this  study. 

REGRESSION  DATA  FOR  THE  INITIAL  LOSS  PARAMETER  (STRTL) 

Because  STRTL  is  storm  dependent,  a  method  of  rela¬ 
ting  the  STRTL  values  to  the  total  amount  of  rain  occurring 
in  the  2,  3,  4,  5,  7,  8,  and  10-day  period  prior  to  the 
storm  was  attempted.  The  statistical  package  MINITAB  was 
used  to  perform  regression  analysis  on  the  data.  The  data 
are  displayed  in  Table  4.3.  The  reader  may  obtain  further 
information  on  MINITAB  by  reading  the  MINITAB  Reference 
Manual  or  by  contacting  MINITAB  PROJECT,  Statistics  Depart¬ 
ment,  215  Pond  Laboratory,  Pennsylvania  State  University, 
University  Park,  PA  16802. 

The  values  in  columns  two  through  eight  were  obtained 
by  summing  the  products  of  the  total  rainfall  recorded 
for  each  rain  gage  covering  a  portion  of  the  basin  and  the 
fraction  of  the  drainage  area  covered  by  that  gage.  The 
analysis  of  this  data  will  be  discussed  in  Chapter  Five. 


I 

I 


Day  Period 


Basin/S  to  rm 

STRTL 

D2 

D3  D4  D5  D7  D8 

DIO 

Butternut 

*  10-18-75 

.31 

.31  .8?  .87  1.25  1-25 

1.25 

10-09-76 

/53 

0 

00000 

.20 

10-21-76 

0 

0 

.01  .01  .40  1.29  1-29 

1.29 

9-26-77 

0 

.73 

.73  -75  1.14  3.19  3-19 

5-40 

10-17-77 

0 

.80 

.8  1.01  1.01  1.73  1-73 

2.02 

Canawacta 

6-26-68 

1.09 

.09 

.09  .09  .16  .38  .38 

.40 

11-18-68 

.26 

.36 

.36  .61  .72  I.47  1.4? 

1.96 

6-24-69 

2.26 

.44 

•51  -51  -51  -81  .81 

.85 

*  7-31-71 

•03 

.88  .88  1.6  1.6  1.6 

2.67 

6-22/23-72 

.08 

0 

0  0  .78  1.67  1.67 

1.67 

7-03-74 

.98 

0 

*  *  FOULED*  * 

* 

Charlotte 

6-29-73 

1.47 

.01 

.01  .14  02  .41  .41 

.5b 

12-21-73 

.38 

.05 

•59  1.54  1.54  1-71  1.71 

1.8 

Corey 

11-14-72 

03 

.46 

.46  .66  .80  .99  .99 

.99 

9-26-75 

.32 

•37 

.46  .46  .62  .77  *77 

.77 

10-09-76 

1.33 

0 

00000 

.06 

5-11-78 

.56 

0 

0  .21  .33  -33  *42 

.87 

*  10-05-79 

.69 

1.14  1.18  1.18  2.02  2.02 

2.02 

11-26-79 

.14 

•  76 

.76  .76  .77  .77  .77 

.77 

6-05-82 

.10 

1.11 

1.11  1.78  1.79  *  * 

• 

Elk  Run 
10-22-70 

1.01 

i  • 

0 

0  0  .09  .55  -55 

1.29 

11-14-72 

•  39 

0 

0  0  1.18  2.28  2.28 

2.28 

9-26-75 

.01 

.45 

.7  .7  .93  1-05  1*05 

1.05 

12-09-76 

1.45 

0 

0  0  0  0  0 

.12 

10-20-76 

.69 

0 

0  0  0  .29  .29 

.29 

5-14-78 

.05 

0 

0  0  .78  .78  .78 

1.55 

*  8-06-78 

2. 

2.  2.  2.4  2.75  2.75 

3.3 

Table  4.3 

Previous  Rainfall  Amounts 
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Day  Period 


Basin/Storm 

STRTL 

D2 

D3 

D4 

D5 

D7 

D8 

DIO 

Five  Mile 

5-07-75 

0 

.1 

* 

Fouled  * 

* 

« 

* 

9-26-75 

•  92 

.6 

.9 

.9 

•  9 

1.1 

1.1 

1.1 

9^14-77 

1.34 

.01 

.01 

.01 

.01 

.01 

.18 

.27 

*  9-25-77 

.01 

.01 

2.09 

2.19 

3-59 

4.6 

504 

Newtown 

9-26-75 

.01 

.81 

H 

00 

• 

-3- 

00 

• 

•  94 

1.43 

1.43 

1.43 

*10-09-76 

•  36 

0 

0 

0 

0 

0 

0 

.1 

10-21-76 

0 

0 

0 

•05 

.26 

.27 

.27 

9-25-77 

.  66 

.04 

.09 

.46 

.80 

1.69 

2.48 

3*9 

11-04-77 

.65 

0 

0 

0 

0 

0 

0 

0 

11-08-77 

•13 

.12 

1.64 

2.05 

2.05 

2.05 

2.05 

2.05 

11-11-77 

.19 

•54 

I.09 

1.09 

1.85 

3.02 

3.02 

3.02 

10-17-77 

.01 

•  25 

•25 

.26 

.27 

.76 

•  76 

•93 

10-06-79 

.96 

•  52 

•52 

.52 

•  52 

1.44 

1.44 

1.44 

11-26-79 

•39 

•3 

.32 

.4 

.4 

.4 

.4 

•  43 

*10-28-81 

.14 

.28 

.28 

.28 

.4 

.56 

•  57 

Ostelic 

7-03-74 

2.33 

.11 

.16 

.47 

•  47 

1.56 

1.56 

1.56 

9-26-75 

1-3 

.24 

.24 

.49 

.49 

.91 

.98 

•  98 

*10-09-76 

0 

0 

0 

0 

0 

0 

•  13 

10-17-77 

.04 

.63 

.63 

.63 

.85 

1.34 

1.82 

2.02 

4-05-78 

0 

.01 

.1 

.26 

.26 

.28 

.48 

1.26 

Owego 

11-08-72 

.61 

.01 

.07 

.08 

•  93 

•  99 

.99 

1.6 

9-26-75 

.29 

.46 

.46 

.63 

.78 

1-15 

1.2 

1.2 

*10-09-76 

.03 

.03 

.03 

.03 

•03 

.03 

•13 

9-25-77 

0 

.03 

•  1 

1.09 

1.45 

2.28 

4.0 

5.25 

10-17-77 

.06 

•3 

•3 

.44 

.45 

1.1 

1.18 

1.37 

Table  4.3.  (Continued) 
Previous  Rainfall  Amounts 
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RESULTS  AND  DISCUSSION 

PURPOSE  AND  SCOPE 

The  purpose  of  this  chapter  is  to  present  the 
results  of  the  calibration  effort  for  the  basins  and 
storms  selected  as  the  control  group;  these  comprise 
selection  of  the  unit  hydrograph  and  regional  loss  rate 
parameters,  determination  of  regression  relationships, 
and  verification  of  the  results  using  uncontaminated 
flood  events .  The  scope  of  this  chapter  is  to  display 
the  results  and  to  discuss  the  findings  in  relation  to 
an  earlier  study. 

MODEL  CALIBRATION  FOR  EACH  BASIN 

The  'best*  fit  values  of  the  four  model  parameters 
for  each  event  are  presented  in  Table  5-1 •  Each  event 
was  initially  optimized  using  the  model’s  default  values. 
Prom  the  initial  results  the  starting  values  of  the 
parameters  were  manipulated  until  the  ’best'  fit  choices 
produced  the  ’best'  statistical  summary  results.  These 
summaries  as  well  as  other  'best'  fit  results  may  be 
found  in  Appendix  C. 

Once  the  optimized  parameters  were  determined  for 
the  control  storms,  the  parameter  for  the  constant  loss 
rate,  CNSTL,  was  averaged  for  each  basin.  For  certain 
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Table  5-1- 

■Best-Fit"  Values  Optimization  Results. 


Basin/Storm 

Butternut  Creek 

10-10-76 
10-21-76 

9-26-77 

10-17-77 


Canasawacta  Creek 
6-26-68 

11-18-68 

6- 24-69 
6-22/23-72 

7- 03-74 

Charlotte  Creek 

6-29-73 

12-21-73 

Corey  Creek 

11-14-72 

9- 26-75 

10- 09-76 
5*14-78 

11- 26-79 
6-05-82 

Elk  Run  Creek 
10-22-70 

11-14-72 
9-26-75 

10- 04-76 
10-20-76 

5-14-78 


TC 

R 

STRTL 

CNSTL 

13-83 

6.24 

.53 

-05 

10.98 

14.06 

0 

0 

9,93 

10.52 

0 

0 

14.77 

4.53 

0 

0 

4.25 

9-69 

1.23 

.06 

2.25 

10-61 

0 

0 

1.83 

12.37 

2.0 

.20 

3-33 

7.94 

-99 

01 

6.13 

4.10 

-9 

.07 

14.61 

20.90 

-92 

-09 

4.21 

15-21 

-50 

.04 

1.62 

10.39 

•30 

.12 

2.27 

2.71 

.39 

.  08 

3.97 

4.73 

1 .62 

-15 

1.03 

6.94 

.28 

.15 

2.60 

3-73 

.22 

.09 

1.05 

11.96 

.04 

.12 

1.07 

5.95 

1.25 

.16 

1.50 

8.54 

.28 

.14 

1 .03 

4.12 

.80 

.10 

4.83 

1.69 

.62 

.07 

1.92 

3-69 

•83 

.07 

1.06 

2  .36 

.08 

.03 

Table  5«1  Continued 
Best-Pit"  Values  Optimization  Results. 


Basin/Storm 


STRTL 


Five  Mile  Creek 


5-07-75 

14.73 

25.44 

.2 

9-26-75 

4.27 

26.76 

.89 

9-14-77 

9.65 

20.35 

1.09 

Newtown  Creek 

9-26-75 

12.48 

20  58 

.03 

10-21-76 

11.01 

4.84 

•  36 

9-25-77 

10.42 

13-59 

.66 

10-17-77 

8.31 

16.02 

.01 

H-04-77 

9.18 

11.90 

•  65 

11-08-77 

8.4l 

14.53 

.13 

11-11-77 

7.86 

13-48 

.19 

10-06-79 

2.64 

18.50 

•  72 

11-26-79 

5.79 

17.17 

.04 

Ostelic  River 

7-03-74 

10.37 

14.26 

2.0 

9-26-75 

6.14 

25.69 

127 

10-17-77 

15-46 

29.53 

.01 

4-05-78 

6.32 

26.52 

.05 

Owego  Creek 

11-08-72 

8.35 

18.94 

.16 

9-26-75 

7.09 

15.23 

.78 

9-25-77 

6.15 

17.78 

.21 

10-17-77 

10.12 

17.28 

.12 
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values  within  the  observed  range  of  this  parameter  (around 
its  average)  and  depending  on  its  performance  the  basin 
value  of  CNSTL  selected  for  each  basin  is  displayed  in 
Table  5*2. 

Because  of  the  unique  structures  of  the  Butternut  and 
Five  Mile  Creeks,  their  values  do  not  appear  to  fit  into 
the  general  pattern  of  increasing  CNSTL  values  as  one 
proceeds  southwestwardly  from  the  Upper  Susquehanna  re¬ 
gion.  Butternut  Creek  has  a  rather  large  value  for  SWT 

(80.26  ft/mi)  and  the  largest  value  for  SF1  (7.27)  implying 
a  long, narrow,  steep-sided  basin.  Such  a  basin  would  not 
allow  as  much  water  to  be  intercepted/infiltrated  as  the 
other  basins.  On  the  other  hand,  a  casual  observation  of 
the  Five  Mile  Creek  basin  on  a  topographic  map  shows 
a  great  deal  of  swamp  area.  This  condition  would  retard 
the  flood  wave;  however,  less  rainfall  would  be  lost 
due  to  the  already  wet  soil  condition. 

Holding  the  parameter  CNSTL  constant  for  each  basin, 
the  other  three  parameters  were  optimized.  The  values  of 
the  initial  loss  parameter,  STRTL,  were  regressed  against 
the  amount  of  total  precipitation  falling  at  different 
time  periods  preceeding  the  storm.  The  best  derived  regres¬ 
sion  equation  was  STRTL  =  0.861  -  0.504  *  D5,  where  D5  is 

the  total  amount  of  rainfall  occurring  five  days  prior  to 
the  beginning  of  the  storm.  See  Figure  5*1  for  a  plot 
of  STRTL  vs  D5  for  all  of  the  control  storms.  This  equa¬ 
tion  has  an  unadjusted  multiple-determination  coefficient, 
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Table  5*2  • 

Basin  Averaged  Constant  Loss  Rates  (CNSTL). 


CNSTL  (inches/hour) 


Butternut  Creek  0.00 
Canasawacta  Creek  0.05 
Charlotte  Creek  0.05 
Corey  Creek  0.10 
Blk  Run  Creek  0.10 
Five  Mile  Creek  0.03 
Newtown  Creek  0.0? 
Ostelic  River  0*05 
Owego  Creek  0.05 
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R^,  equal  to  .192. 


2.40+ 


** 


1.80+ 


-  * 
—  2 
1.20+ 


* 


* 


* 

* 


Although  the  coefficient  is  not  as  high  as  desired, 
the  D5  term  gave  the  best  correlation.  The  control  storms 
plotted  below  the  line  in  Figure  5*1  will  give  conserva¬ 
tive  forecasts.  This  means  higher  computed  peaks  than 
the  observed  peaks.  The  STRTL  terra  is  rather  sensitive 
(see  Sensitivity  Analysis  section,  later),  however,  in 
the  absence  of  any  other  method, this  approach  was  used. 

The  values  for  STRTL  were  predicted  by  means  of  this 
regression  equation.  Each  storm  was  then  optimized  for 
the  TC  and  R  parameters .  The  final  parameter  values  for 
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the  control  storms  are  listed  in  Table  5 >3* 

These  values  were  used  in  computing  and  comparing 
the  verification  storm  group  hydrographs  with  the  observed 
flood  hydrographs.  The  statistical  results  are  presented 
in  Appendix  D  for  the  verification  storms. 

COMPARISON  WITH  PARAMETERS  DETERMINED  PRIOR  TO  HURRICANE 

AGNES  (1972) 

In  June  1970  the  Susquehanna  River  Basin  Study  Coor¬ 
dinating  Committee  published  an  extensive  work  on  the 
Susquehanna  basin.  Appendix  D  (Hydrology)  of  that  study 
contained  the  Clark  unit  hydrograph  parameters  for  several 
gaged  basins,  obtained  on  the  basis  of  a  similar  computer 
program  developed  by  the  Hydrologic  Engineering  Center, 

U  S.  Army  Engineer  District,  Sacramento,  CA,  (now  HEC  is 
located  in  Davis,  CA). 

The  data  base  used  in  the  study  consisted  of  four 
major  flood  events.  Basin  average  rainfall  was  determined 
by  Thiessen  polygons .  Two  or  four-hour  unit  hydrographs 
were  developed  depending  on  the  drainage,  area  of  the  basin. 
The  TC  and  R  parameter  values  are  presented  in  Table  5 -4 • 
Not  all  of  the  basins  could  be  compared  with  the  current 
results  because  some  stream  gages  have  been  relocated 
since  the  publication  of  the  study.  Table  5.4  shows  a 

very  good  comparison  between  most  of  the  terms.  The  values 
of  R  for  Charlotte  and  Five  Mile  Creeks  are  not  in  good 
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Table  5*3* 
Derived  Parameters. 


Basin 

TC 

R  (hours) 

Butternut  Creek 

11.66 

6.76 

Canasawacta  Creek 

8.65 

8.59 

Charlotte  Creek 

11.68 

27-70 

Corey  Creek 

1.90 

8.27 

Elk  Run  Creek 

2.02 

8.13 

Five  Mile  Creek 

7.47 

46.12 

Newtown  Creek 

9.00 

14.46 

Ostelic  River 

9.80 

36.13 

Owego  Creek 

9.35 

13.83 
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agreement.  A  possible  explanation  for  the  Charlotte  Creek 
discrepancy  is  the  fact  that  only  two  storms  meeting  the 
established  criteria  in  the  past  fifteen  years  could  be 
found  for  this  basin.  The  expected  value  variance  of  any 
derived  parameter  based  on  only  two  data  points  is  very 
great.  Therefore,  one  would  not  expect  the  values  of  the 
two  works  to  be  very  close. 

The  discrepancy  of  the  R  values  for  Five  Mile  Creek 
could  be  related  to  the  goodness-of-fit  of  the  storms 
selected  for  this  study.  Looking  at  the  plots  of  the 
'best*  fit  computed  flow  with  the  observed  flow,  the  reader 
will  notice  that  the  computed  hydrograph  had  a  rather 
difficult  time  matching  the  observed  hydrograph.  Whether 
the  rain  gage  in  the  middle  of  the  basin  failed  to  provide 
accurate  data  or  the  U.S.G.S.  gaging  station  measurements 
were  fouled  by  mechanical  defects  is  not  clear.  The  author 
is  inclined  to  doubt  the  observed  hydrograph  measurement 
because  of  the  peculiar  combination  of  spike  and  round 
portions . 

Study (Jun  70)  Present  Results 


TC 

R 

TC 

R 

Canasawacta  Ck. 

2.03 

11.04 

8.65 

8.59 

Charlotte  Ck. 

II.83 

14.79 

11.68 

27.70 

Five  Mile  Creek 

5.50 

19-00 

7.47 

46.12 

Newtown  Creek 

971 

18.39 

9.00 

14.46 

Owego  Creek 

10.10 

11.82 

905 

13-83 

Table  5-^. 

Parameter  Comparison  with  Previous  Study 


SENSITIVITY  ANALYSIS  OP  PARAMETERS 
A  sensitivity  analysis  was  performed  on  a  selected 
storm  whose  ’best'  fit  results  initially  yielded  a  very 
good  statistical  summary.  The  storm  selected  from  the  con¬ 
trol  group  was  the  Newtown  Creek  11-04-77  storm.  Holding 
the  other  three  parameters  constant,  the  parameter  TC 
was  decreased  and  increased  5  and  10 %.  The  same  was  done 
for  the  parameter  R.  The  loss  rate  parameters  are  not 
normally  known  to  such  a  high  degree  of  accuracy.  There¬ 
fore,  the  loss  rate  parameters  were  not  varied  by  5  and 
10%  because  the  results  were  not  significantly  different. 
For  the  initial  loss  rate  parameter,  STRTL,  the  value 
was  decreased  and  increased  .1  and  .2  incnes  to  show  a 
better  variability.  The  parameter  for  constant  loss  rate, 
CNSTL,  was  varied  by  .02  and  .04  inches/hour  each  way. 

These  variations  are  reasonable  for  any  basin  under  inves¬ 
tigation.  The  results  are  displayed  in  Table  5*5* 

The  reader  is  directed  to  the  Peak  Flow  column.  The 
Peak  Flow  percent  was  not  very  sensitive  to  variations 
in  the  TC  and  R  terms.  The  most  change  was  about  6%  between 
the  computed  and  observed  peaks.  Variations  in  the  normal 
range  of  expected  values  of  the  loss  parameters  are  a  dif¬ 
ferent  story.  A  decrease  of  STRTL  of  .2  inch  caused  an 
overestimation  of  the  peak  by  lk%.  While  a  .04  inciv/hour 
decrease  or  increase  in  the  CNSTL  produced  approximately 
a  40Jt  change  in  both  directions. 


Sum  of  Equiv  Mean  Time  To  Lag  #  Time  AVG 
Flows  Depth  Flow  C.M.  C.M.  to  Peak  of  %  ABS 
Adjustments  Diff.Diff .  Diff  Diff  C.M.  *  Flow  Peak  Error  OF 
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Sensitivity  Results 


The  results  of  the  analysis  show  that  the  loss  parame¬ 
ters  are  the  most  sensitive  and  must  be  estimated  accurately 
to  give  reasonable  results . 

The  reader  may  recall  the  discussion  in  Chapter  Three 
concerning  the  optimized  solutions  of  HEC-1  not  being  the 
global  solution.  This  point  is  reflected  in  this  analysis 
where  a  5#  increase  in  TC  produced  a  better  fit  than  the 
original  controlled  parameter  storm. 

REGRESSION  ANALYSIS  OF  THE  PARAMETERS  TC  AND  R 
Various  combinations  of  the  basin  characteristics 
listed  in  Table  4.  were  tried.  The  only  restriction  placed 
on  the  combinations  was  that  only  one  slope  value  and 
only  one  shape  factor  were  permitted  in  any  one  combina¬ 
tion.  Combinations  of  the  characteristics  as  well  as 
their  log  values  were  evaluated.  Table  5.6  lists  the 
top  combinations  of  each  type  along  with  the  statistical 
coefficients.  Table  5-7  lists  some  of  the  characteristic 
combinations  used. 


7.942.63(DA)*66/(ALT)1’44 

rL 

r2-ADJ 

Standard 

Deviation 

TC  * 

.910 

.879 

0.2474 

TC  « 

450. 34/(S2) 1,08 

.878 

.861 

0.2659 

TC  « 

5-93  (DA/S)-412 

.871 

.853 

0.2736 

TC  « 

1.73(tVDA/S) 

.861 

.842 

0.2836 

TC  * 

3.60  (I/^V716 

5* 

CO 

e 

.825 

0.2980 

Table  5^ •  Highest  Correlated  Regression  Equations. 
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0 

O 

Standard 

sf 

R-ADJ 

Deviation 

TC+R  =626.41/(S)*87 

.731 

•693 

0.3365 

TC+R=l4.15(L/K3)  *549 

.686 

•  64l 

O.3638 

TC+R=20 . 9 1 ( DA/S ) * 31 2 

.686 

.641 

O.3639 

TC+R=8 . 25 ( L  V&A/S ) * 35 

.677 

.630 

O.3693 

TC+R=15.18(L/\§WT) *653 

.667 

.619 

0.3748 

TC+R=507 . ?6/(S2 ) * 794 

.658 

.610 

0.3795 

Table  5  >6  . 

Highest  Correlated  Regression  Equations 

(Continued) 

DA/S 

DA/SWT  L  LCM/VS  lVda/S 

L/VSf  L/VSWT 

Butternut  Ck .  1.75 

•  74 

30.53  27.55 

3-57 

2.33 

Canasawacta  Ck .  1.24 

.84 

15.44  13.85 

1.82 

1.50 

Charlotte  Ck.  4.83 

2.4? 

63.80  54.42 

4.21 

3.01 

Corey  Ck.  .09 

.11 

1.4  I.61 

.46 

.52 

Elk  Run  Ck.  .10 

.08 

2.91  2.10 

.66 

.59 

Five  Mile  Ck.  2.25 

1.00 

37.27  29.23 

3-58 

2.38 

Newtown  Ck .  1.72 

1.10 

28.48  23.48 

2.67 

2.13 

Ostelic  Ck.  6.2 

2.84 

17.72  73.46 

6.06 

4.10 

Owego  Ck .  6.33 

3.40 

90.20  76.92 

5.66 

4.15 

Table  5-7  • 

Basin  Parameter  Combinations 
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The  first  equation  in  each  group  in  the  table  has  the 
highest  correlation  coefficients,  lowest  variance,  and  least 
number  of  parameters.  The  parameters  relating  best  to  the  Clark 
unit  hydrograph  parameters  TC  and  R  were  found  to  be  the  drainage 
area  (DA),  altitude  difference  (ALT),  and  the  main  channel  slope 
(S).  The  shape  factors  were  found  to  be  highly  related  to  the 
other  terms.  They  were  eliminated  without  exception  in  every 
stepwise  regression  operation.  The  term  S2  faired  better  in  the 
TC  group  than  the  S  term.  The  opposite  is  true  in  the  TC  +  R 
group.  However,  the  difference  between  them  is  not  great  and 
does  not  warrant  the  additional  effort  of  deriving  S2. 

AUTOMATIC  CALIBRATION  RESULTS 
In  Chapter  Three  under  the  section  "Calibration  of  Model 
Parameters",  it  was  mentioned  that  the  optimization  technique 
employed  by  HEC-1  does  not  guarantee  that  the  parameters  derived 
are  global  solutions.  An  example  of  this  is  displayed  in  Figures 
5.2  and  5.3.  Figure  5*2  is  the  graphical  results  obtained  using 
the  HEC-1  initial  parameters  as  starting  points.  Figures  5*3 
is  the  final  results  used  in  this  analysis  for  the  Butternut 
Creek  (9-26-77)  storm. 
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Figure  5*3*  Finial  Optimization  Results 
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CHAPTER  SIX 


SUMMARY  AND  CONCLUSIONS 


McSparran  (1968,  p.937)  found  in  his  study  of  32 
watersheds  varying  in  drainage  area  from  2  to  210  square 
miles  in  Pennsylvania  that  relating  hydrograph  parameters 
to  watershed  characteristics  is  not  an  easy  task.  He  found 
nonlinearity  factors  present  in  every  storm,  adjusting 
and  altering  the  parameters.  As  a  result  he  felt  that 
the  unit  hydrograph  method  was  not  as  accurate  as  it 
should  be.  He  stated  that  major  floods  would  be  underpre¬ 
dicted  and  on  the  average  the  observed  peaks  would  be 
underpredicted  by  22#  using  the  Snyder  unit  hydrograph 
method.  Yet  by  using  the  Soil  Conservation  Service  (SCS) 
method  the  average  peaks  were  predicted  51#  too  high. 

On  the  other  hand,  Heerdegen  (1974,  p.ll43)  concluded 
with  a  resounding  "affirmative"  to  his  article's  question- 
title,  "Unit  Hydrograph i  A  Satisfactory  Model  of  Watershed 
Response?"  Sixty  measured  parameters  for  each  of  the  six¬ 
teen  watersheds  in  Pennsylvania  were  used  in  analyzing 
ninety-six  flood  events.  Five  parameters  were  eventually 
selected  as  the  most  influential  for  his  area.  However, 
they  were  not  constant  from  area  to  area.  In  his  study 
he  concluded  that  the  storm  variable  of  precipitation 
duration  and  volume  of  precipitation  excess  are  major 
factors  in  flood  hydrograph  responses. 
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SOURCES  OP  ERROR 

Whenever  one  attempts  to  measure  a  physical  property 
in  nature,  there  is  some  error.  Rain  gages  and  streamflow 
gages  leave  a  lot  to  be  desired  when  it  comes  to  accurate 
measurements.  Rain  gages  are  point  measuring  instruments 
and  they  hardly  ever  represent  the  amount  of  rainfall 
deposited  on  a  basin  accurately.  In  some  cases  rain 
gages  at  a  great  distance  from  the  basins  had  to  be  used, 
because  the  local  data  were  fouled  or  not  reported. 

The  choice  of  model  is  a  source  of  error.  Every 
model  makes  some  simplifying  assumptions  for  ease  of 
computation,  cost,  or  speed.  HEC-1  is  a  lumped  or  bulk 
parameter  model.  It  groups  all  of  the  irregularities 
and  unique  features  of  a  basin  into  one  big  box.  Richards 
(1955.  p-110)  presented  nine  factors  affecting  runoff. 

Some  of  these  factors,  such  as  moving  storms  and  intensity 
variations,  are  ignored  by  HEC-1. 

In  the  absence  of  specific  information,  a  synthetic 
time  4rea  curve  (supplied  by  HEC-1)  was  made  to  create  a 
translation  hydrograph  for  each  basin.  The  basins  varied 
in  many  sizes  and  shapes  and  forcing  the  time-area  curve 
to  fit  all  the  basins  seems  unrealistic. 

The  type  of  optimization  routine,  i.e.  univariant 
search  technique,  does  not  guarantee  the  results  to  be 
a  global  solution  and  indeed  the  user  of  the  model  must 
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be  aware  of  its  hazards.  Picking  new  starting  points  and 
adjusting  the  parameters  to  get  a  better  fit  is  sometimes 
more  luck  than  skill. 

The  loss  rate  selection  is  a  major  source  of  error. 

The  inability  to  obtain  a  consistent  and  reliable  STRTL 
predicting  device  greatly  influenced  the  final  results 
and  future  applications  of  this  study.  Bulk  loss  rates 
are  perhaps  better  suited  for  large  areas. 

Selection  of  the  Clark  method  to  derive  flood  hydro¬ 
graphs  affected  the  results.  The  application  of  a  linear 
numerical  model,  although  frequently  used  in  government 
and  private  studies,  to  a  nonlinear  problem  induces  error. 

Are  the  areas  of  interest,  i.e.  the  Upper  Susquehanna 
and  Chemung  River  Basins,  compatable?  The  nine  stream 
basins  are  supposed  to  be  hydrologically  similar.  Are 
they?  What  does  hydrologically  similar  mean?  This  author 
concluded  that  the  basins  are  fairly  hydrologically  similar 
because  the  slopes  and  valley  structures  appear  to  be 
consistent  among  the  nine  basins  and  most  major  storms 
track  across  both  river  basins.  What  basin  or  climatic 
factors  were  not  included  which  would  have  improved  the 
correlation  between  the  model  parameters  and  the  basin 
characteristics?  Chow  (1964,  p. 27-104)  states  that  some 
areas  are  just  not  adequate  for  unit  hydrographs. 

The  unit  hydrograph  assumptions  of  time  invariance 
and  linear  runoff  responses  are  not  strictly  valid  in 
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real  life.  Basin  characteristics  do  change  with  time,  e.g. 
land  cover  and  seasonal  changes  influence  loss  rates. 

The  selection  of  storms  will  affect  the  results. 
Although  the  storms  are  the  largest  occurring  in  the  last 
ten  to  fifteen  years  (unaffected  by  snow),  most  are  not 
to  be  considered  major  storms.  Sokolov  (1976,  p.199)  quoted 
U.S.  Army  Corps  of  Engineers'  studies  when  they  said  that 
in  the  majority  of  the  basins  studied  the  peak  ordinates 
of  the  unit  hydrographs  derived  from  a  major  flood  are 
generally  significantly  greater  than  those  derived  from 
data  from  a  minor  flood.  Those  studies  stated  that  unit 
hydrograph  peak  ordinates  from  major  floods  were  twenty- 

five  to  fifty  percent  greater  than  those  derived  from  minor 
floods.  Therefore,  one  would  expect  the  results  obtained 
in  the  application  of  the  regression  relationships  to  be 
somewhat  underpredictive,  i.e.  on  the  conservative  side. 

The  storms  selected  cover  several  seasons  from  early 
April  until  December.  Although  Heerdegen  (1968,  p.ll43) 
found  no  significant  seasonal  differentiation  in  the  unit 
hydrograph  at  the  five  percent  level,  Sherman  (1932,  p-504) 
felt  that  unit  hydrographs  should  be  seasonally  segregated. 

LIMITS  OF  APPLICATION 

The  data  were  obtained  from  streams  and  rain  gages 
in  south  central  New  York  and  north  central  Pennsylvania 
and  they  should  be  applied  in  that  area .  With  the  exception 
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of  Five  Mile  Creek,  the  streams  appear  to  be  free  from 
a  great  deal  of  valley  storage  and  swampy  areas.  The  results 
were  obtained  without  snowmelt  or  ice  jams  affecting  the 
outlet  flow  rate.  Therefore,  the  results  should  apply 
primarily  to  summer  and  fall  seasons. 

AREAS  IN  NEED  OF  FUTURE  INVESTIGATION 
The  lack  of  adequate  soil  moisture  data  hindered 
the  proper  calibration  of  the  parameter  STRTL.  Without 
adequate  or  reasonable  loss  rate  parameter  values,  the 
model  will  not  serve  as  a  reliable  predictive  tool. 


SUMMARY 

By  using  the  Flood  Hydrograph,  HEC-1,  computer 
package,  the  Clark  unit  hydrograph,  and  the  initial/ 
uniform  loss  rate,  parameters  were  determined  for  nine 
watersheds  in  the  Chemung  and  Upper  Susquehanna  River 
basins.  Fifty- two  storms  were  used  in  the  analysis  on 
the  watersheds  which  range  in  drainage  areas  from  10.2 
to  185  square  miles.  Regression  equations  predicting 
the  time  of  concentration,  TC,  and  the  term  TC  +  R 
were  derived. 
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CONCLUSIONS 

Specifically,  the  conclusions  of  this  thesis  arei 

1.  HEC-1  may  be  applied  successfully  in  the  study  area. 

2.  Sensitivity  of  the  parameters  TC,  R,  STRTL,  and 
CNSTL  in  HEC-1  are  established. 

3.  The  HEC-1  optimization  search  technique  should 
be  improved  to  give  a  global  solution. 

b.  The  relationship  to  determine  the  term  STRTL  needs 
to  be  improved. 

5.  Regional  values  of  the  term  CNSTL  are  now  established 

6.  The  initial  and  constant  loss  rate  method  is  not 
the  best  method  but  it  works  fairly  well. 

7.  The  results  should  be  applied  seasonally. 

8.  Good  correlation  between  the  parameters  TC  and 
R  with  basin  characteristics  were  obtained.  The 
recommended  equations  are  the  first  in  each  group 
in  Figure  5* 

9*  A  very  good  comparison  with  a  1970  study  was 
obtained  for  the  terms  TC  and  R. 


Appendix  A  Watershed  Rain  Gage  Locations  snd  Thiessen 
Polygons. 

Recording  Stations!  Edmeston 

Edmeston  3N 

Non-recording  Stations!  Cooperstown 

New  Berlin 


i  / 


[«-•  with  Edmeston  3N 
— without 
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Recording  Stations  i  Plymouth 

Plymouth  2SSE 

Non-recording  Stations*  Cincinatus 

i  Lincklaen 
Norwich 
Norwich  1NE 
Sherburne 


CANASAWACTA  CREEK  near 
Plymouth,  NY 


Cobleskill  2 


Oneonta 

State 

Univ. 


Recording  Stations i  Canton 

Jackson  Summit 

Non-recording  Stations i  Canton 

Covington 

Troy 


•  Jackson  Summit 


EIK  RUN  CREEK 


COREY  CREEK  near 
Mains burg,  PA 


Canton 


Recording  Stations i  None 

Non-recording  Stations t  Bath 

Pratt burg 


Recording  Stations! 


NEWTOWN  CREEK  at 
Elmira,  NY 


Recording  Stations i  Plymouth 

Plymouth  2S3E 


Non- recording  Sta.si  Cincinatus 

DeRuyter 

Lincklaen 


OTSELIC  RIVER  at 
Cincinatus,  I'Tf 
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•  Cortland 


Recording  Stations t 
Hunt’s  Comers 
Mt.  Pleasent  Farm 
Newark  Valley 


Non-recording  Stations « 
Candor 
Cortland 
Lisle 


Appendix  B 

Input  Card  Deck  Description 

A  brief  line-by-line  explanation  of  a  typical  card 
deck  used  in  this  study  follows i 

ID  —  (Identification  card)  Serves  to  identify  job  by 
providing  an  opportunity  to  make  comments  or  to  describe  an 
operation  to  be  performed. 

IT  —  (Identify  Time  periods  card)  Serves  to  establish 
hydrograph  time  axis  (60  minutes),  start  date  of  storm 
hydrograph  (24  September  19??),  start  time  of  storm 
hydrograph  (1300  hours),  and  number  of  hydrograph  or¬ 
dinates  which  must  match  the  total  number  of  flowrate 
entries  on  the  QO  cards. 

10  —  (Identify  Output  formats)  Serves  to  request  supression 
or  printing  of  selected  output,  calculations,  or  plots. 

OU  —  (Optimization  limits  for  Unit  hydrographs)  The  first 
number  tells  the  program  at  which  flood  hydrograph  ordinates 
optimization  must  be  started.  The  second  establishes  the 
optimization  ending  ordinate. 

PG  —  (Precipitation  Gage)  Serves  to  identify  rainfall 
gages.  If  a  number  follows  the  gage  name,  the  gage  is 
non-recording  and  the  number  is  the  total  storm  rainfall 
in  inches.  If  no  amount  of  rainfall  is  given,  then  the 
gage  is  a  recording  gage  and  must  be  followed  by  a  PI  card. 


AD-A129  964  WATERSHED  MODELLING  IN  THE  CHEMUNG  AND  UPPER 

SUSQUEHANNA  RIVER  BASINS  USI..(U)  ARMY  MILITARY 
PERSONNEL  CENTER  ALEXANDRIA  VA  J  C  STYRON  20  MAY  83 

F/G  8/8 


\ 


UNCLASSIFIED 


NL 


Ill  1.0  SB*  Itt 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL,  BUREAU  OF  STANDARDS- I963-A 


8U 


PI  —  (Precipitation,  Incremental)  Serves  to  follow  a 
recording  rain  gage  PG  card  to  describe  the  temporal 
rainfall  distribution  of  the  recording  gage.  The  values 
are  spaced  at  the  time  interval  established  on  the  IT 
card.  The  first  value  must  represent  the  amount  of  rain 
falling  in  the  first  time  period  indicated  on  the  IT  card. 
The  start  of  the  rain  must  coincide  with  the  start  of  the 
flood  hydrograph  flowrate  on  the  QO  cards. 

KK  —  (Kontrol  Kard)  Serves  to  identify  the  name  of  the 
stream  gaging  station  at  the  basin  outlet  and  to  separate 
the  precipitation  data  cards  from  the  outflow  data. 

00  —  (Discharge  ( Q )  Ordinate)  Serves  to  input  the  observed 
hydrograph  ordinates.  The  first  flow  value  must  match  the 
first  PI  card  entry  as  well  as  the  IT  card,  i.e.  all 
occur  at  the  same  instant. 

PT  —  (Precipitation  Gage  Total)  Serves  to  identify  those 
precipitation  gages  (either  recording  or  nonrecording) 
which  cover  some  fraction  of  the  basin  and  contribute 
to  the  basin  runoff. 

PW  —  (Precipitation  Gage  Weightings)  Serves  to  provide 
the  fraction  of  the  basin  covered  by  the  gage  above  it. 

PR  —  (Precipitation  Recording)  Serves  to  identify  the 
recording  stations  which  are  used  in  deriving  the  temporal 
rainfall  distribution  pattern. 

PW  —  Same  as  the  PW  card  above  except  it  represents  the 
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fraction  of  the  basin  covered  only  by  recording  stations. 
BA  —  (Basin  Area)  The  first  value  is  the  drainage  area 
while  the  second  is  the  mean  annual  rainfall  for  the 
basin,  if  known. 

BP  —  (Base  Plow)  The  values  of  STRTQ,  QRCSN,  and  RTIOR 
respectively. 

UC  —  (Unit  Hydrograph,  Clark  Method)  The  two  parameters 
are  the  time  of  concentration  (TC)  and  storage  (R), 
respectively.  If  the  numbers  are  positive,  they  remain 
constant  throughout  the  calculations.  If  negative,  the 
optimization  routine  begins  with  these  values  and  then 
adjusts  their  values  after  removing  the  negative  sign. 

If  one  or  both  are  -1,  then  the  optimization  routine 
supplies  an  initial  default  value  and  then  optimizes. 

In  this  model  both  parameters  must  be  either  positive 
or  negative  because  of  their  interrelationship. 

UJ  —  (Loss  Rate,  Initial  and  Uniform)  The  two  parameters 
are  the  initial  loss  (STRTL)  and  the  constant  loss  rate 
(CNSTL),  respectively.  The  same  definitions  for  positive 
and  negative  values  apply  as  above,  except  that  both 
do  not  have  to  be  positive  or  negative  at  the  same  time. 
ZZ  —  (End  of  Data)  HEC-1  begins  calculations  after  all 
data  is  read  into  its  scratch  files. 
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Appendix  D  —  Statistical  Results  of  Verification  Storms 


Basin/Storm 

C.M. 
to  C  .M . 

$6  Peak 
Flow 

Time 
to  Peak 

Avg  % 

Abs  Error 

OBJ 

Func 

Butternut  Ck. 
10-18-75 

1.37 

-14.92 

0. 

15.39 

248 

Canasawacta  Ck. 
7-31-71 

5-73 

130.85 

3. 

238.17 

1514 

Corey  Ck. 
10-09-76 

3.35 

9.79 

-2. 

37.96 

62 

Elk  Run  Ck. 
8-06-78 

44.73 

109.13 

0. 

118.80 

197 

Five  Mile  Ck. 
9-25-77 

-10.92 

-40.37 

2. 

31-17 

529 

Newtown  Ck. 
10-09-76 

1.17 

28.91 

-1. 

32.62 

450 

10-28-81 

38.41 

9.49 

0. 

37.62 

331 

Ostelic  River 
10-09-76 

9.85 

-20.67 

5- 

27.60 

1223 

Owego  Ck. 
10-09-76 

3.15 

36.59 

0. 

47.24 

1903 

29k 


BIBLIOGRAPHY 


Benjamin*  J.R.  and  Cornell,  C. A. (1970). Probability,  Statis¬ 
tics  and  Decision  for  Civil  Engineers.  McGraw-Hill, 

New  York,  419  pp. 

Biswas,  A. K.( 1967). Hydrologic  Engineering  Prior  to  600  B.C. 
Proceedings  ASCE.  Journal  of  Hydraulic  Division.  93. 125. 

Brater,  E.F.( 1940) .The  Unit  Hydrograph  Principle  Appied  to 
Small  Watersheds.  Trans  ASCE.  105 . 1154  pp. 

Chow, V .T . (1964) .  Handbook  of  Applied  Hydrology.  McGraw-Hill, 
New  York,  1 4- 25pp. 

Clark, C .0 . (1945) .Storage  and  the  Unit  Hydrograph.  Amer.  Soc. 
Civ.  Engin.,110,l4l6-l446. 

Cruff,  R.W.and  Rantz, S .E . (1965) •  A  Comparison  of  Methods 
Used  in  Flood- Frequency  Studies  for  Coastal  Basins  in 
California.  Flood  Hydrology,  USGS  Water  Supply  Paper 
1580. 

Dooge,  J.C .1  .(1959) •  A  General  Theory  of  the  Unit  Hydrograph. 
Journal  of  Geophysical  Research.  64,2,241-256. 

Dooge,  J.C. I. (1973).  The  Linear  Theory  of  Hydrologic  Systems. 
U.S.Dept.  of  Agriculture  Tech.  Bulletin  No.  1468.  U .S . 
Govt.  Printing  Office. 

Heerdegen, R.G . (1974) .The  Unit  Hydrograph  1  A  Satisfactory 
Model  of  Watershed  Response?  Water  Resources  Bulletin. 
Amer.  Water  Resources  Assoc . ,10,6, 1143-1159. 

Hoyt,  W.G.  and  Lang be in. W.B. (1955) .Floods. Princeton  Univer¬ 
sity  Press, lip. 

Kinnison,  H.B.and  Colby, B .R . (1945 ) .  Flood  Formulas  Based  on 
Drainage  Basin  Characteristics. Trans  ASCE. 110. 

Langbein,  W.B.,  et  al.  (1947).  Topographic  Characteristics 
of  Drainage  Basins.  USGS  Water  Supply  Paper  968-C. 
125,131.138. 

Linsley,  R,K.  Jr.,  and  Kohler,  M.A.,  and  Paulhus,  J.L.H. 
(1982).  Hydrology  for  Engineers.  Mc-Graw-Hill,  New 
*ork,  339-353  PP- 


295 


296 


McSparran,  J.E.  (1968).  Design  of  Hydrograph  for  Pennsyl¬ 
vania  Watersheds.  Journal  of  the  Hydraulics  Division. 
Proceedings  ASCE,  937-959pp. 

Morgan,  R.,  and  Hullinhors,  D.W.  (1939).  Unit  Hydrographs 

for  Gaged  and  Ungaged  Watersheds,  U.S.  Engineers  Office, 
Binghamton,  New  York. 

Nash,  J.E.  (1957)*  The  Form  of  the  Instantaneous  Hydrograph. 
Intern  Assoc  Sci  Hydrology.  Pub  45,  Vol.  3,  pp. 114-121. 

O'Kelley,  J.J.(l955)*  The  Employment  of  Unit-Hydrographs  to 
Determine  the  Flows  of  Irish  Arterial  Drainage  Channels. 
Proc  Inst  Civil  Engineers,  4,pt.3, 365-4l2pp. 

Richards,  B.D. (1955). Flood  Estimation  and  Control.  Chapman 
and  Hall,  Ltd,  London,  3rd  Ed,  l-8pp. 

Rodriguez- Iturbe,  I.,  Sanabria,  M.G.  and  Caamano,  G.(1982). 
On  the  Climatic  Dependence  of  the  IUHi  A  Rainfall — 
Runoff  Analysis  of  the  Nash  Model  and  the  Geomorpho- 
climatic  Theory.  Water  Resour.  Res. ,18,4.887-903. 

Rodriguez -Iturbe, I .  and  Gonzalez-Sanabria ,  M.(1982).  A  Geo- 
mo  rpho climatic  Theory  of  the  Instantaneous  Unit  Hydro- 
graph.  Water  Resour.  Res, 18, 4  ,877-886  pp. 


Sherman,  L.K.(1932).  Stream  Flow  from  Rainfall  by  the1  . 
unit-graph  method.  Engineering  News  Record.  108, 
501-505. 

Sokolov,  A. A.,  Rantz,  S.E.  and  Roche,  M.(1976).  Flood  flow 
Computation — Methods  Compiled  from  World  Experience, 
The  UNESCO  Press,  171pp. 

Snyder,  F.F.  (1938). Synthetic  Unit  Graphs,  Trans.  Amer. 
Geophysical  Union.  19.  p.447-454. 

Taylor,  A.B.  and  Schwarz,  H.E.(1952).  Unit  Hydrograph  Lay 
and  Peak  Flow  Related  to  Basin  Characteristics,  Trans 
Am  Geophysical  Union. 33.  235-246. 

U.S. Army  Corps  of  Engineers. ( 1981 ) .HEC-1  Flood  Hydrograph 
Package,  Users  Manual,  p.1-80. 

U.S. Army  Corps  of  Engineers. (1982) .Hydrologic  Analysis  of 
Ungaged  Watersheds  Using  HEC-l,  Training  Document 
No.  15,  p. 1-105. 


